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We present two-dimensional (2D) two-particle angular correlations on relative pseudorapidity 
rj and azimuth <j) for charged particles from Au-Au collisions at ^/snn = 62 and 200 GeV with 
transverse momentum pt > 0.15 GeV/c, \r)\ < 1 and 2n azimuth. Observed correlations include 
a same-side (relative azimuth < n/2) 2D peak, a closely-related away-side azimuth dipole, and 
an azimuth quadrupole conventionally associated with elliptic flow. The same-side 2D peak and 
away-side dipole are explained by semihard parton scattering and fragmentation (minijets) in proton- 
proton and peripheral nucleus-nucleus collisions. Those structures follow N-N binary-collision scaling 
in Au-Au collisions until mid-centrality, where a transition to a qualitatively different centrality 
trend occurs within a small centrality interval. Above the transition point the number of same-side 
and away-side correlated pairs increases rapidly relative to binary-collision scaling, the rj width of 
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the same-side 2D peak also increases rapidly (n elongation) and the <j> width actually decreases 
significantly. Those centrality trends are in marked contrast with conventional expectations for 
jet quenching in a dense medium. The observed centrality trends are compared to perturbative 
QCD predictions computed in hijing, which serve as a theoretical baseline, and to the expected 
trends for semihard parton scattering and fragmentation in a thermalized opaque medium predicted 
by theoretical calculations and phenomenological models. We are unable to reconcile a semihard 
parton scattering and fragmentation origin for the observed correlation structure and centrality 
trends with heavy ion collision scenarios which invoke rapid parton thermalization. If the collision 
system turns out to be effectively opaque to few-GeV partons the present observations would be 
inconsistent with the minijet picture discussed here. 

PACS numbers: 25.75.-q, 25.75.Gz 



I. INTRODUCTION 

Many conventional theory descriptions of central colli- 
sions at the Relativistic Heavy Ion Collider (RHIC) full 

5 energy invoke the basic assumption that copious parton 
(mainly gluon) production during initial nucleus-nucleus 
(A-A) contact and subsequent parton rescattering lead 
to a color-deconfined, locally-thermalized quark-gluon 
plasma [l|, 0]. Hydrodynamic models claims of 

10 "perfect liquid" formation 0-QIij], and the relevance of 55 
lattice QCD predictions to RHIC data all rely on as- 
sumed formation of a rapidly-thermalized QCD medium. 
However, experimental confirmation of that assumption 
remains an open question. Although the constituents 

15 of the system may interact strongly, thermalized mat- 60 
ter may not emerge in the time available in relativistic 
collisions . Experimental study of possible rapid ther- 
malization is one of the goals of this paper. 

RHIC heavy ion collisions are studied as a function of 65 
20 nucleus size A, collision energy and centrality to search 
for evidence that an approximately linear superposition 
of nucleon-nucleon (N-N) interactions [l2| expected for 
peripheral A-A collisions evolves with increasing size, 
energy and centrality to a collective system of dense, 70 
25 strongly-interacting QCD matter. In reports by the 
four RHIC experiments [13l4l6l | it was argued that ob- 
servations are consistent with a collective thermalized 
medium. 

High-p t jet tomography was proposed to probe the 75 

30 conjectured QCD medium. Hard-scattered partons pro- 
duced in large-Q interactions during initial A-A con- 
tact [where Q is the parton (actually dijet) energy scale] 
are nominally well-understood probes of collision dy- 
namics and QCD medium pro pert ies (i.e., described by so 

35 perturbative QCD or pQCD) [17j . The underlying as- 
sumption is that formation of a QCD medium should 
modify parton scattering and fragmentation to hadrons 
and may thereby produce deviations of corresponding 
hadron distributions (single-particle spectra and correla- 85 

40 tions) from binary-collision scaling fl3l . [lij . Much atten- 
tion has therefore been paid to high-pt systematics (e.g., 
reduced high-p t hadron yields suppression of jet- 

rclatcd away-side azimuth correlations [19() interpreted 
to reveal strong parton energy loss [l?} • But those results 90 

45 do not distinguish thermalization scenarios from other 



possibilities [IT| . 

In this paper we utilize two-particle angular correla- 
tions among all accepted charged particles and focus on 
those structures associated with semihard parton scat- 
tering and fragmentation [2fjj |. referred to as minijet an- 
gular correlations. Those structures provide a comple- 
mentary approach to medium studies. Inference of jet 
structure (minijcts) from minimum-bias (all particles in 
the pt acceptance) angular correlations [2ll - [24| differs 
qualitatively from high-p t jet methods in that the mini- 
jet analysis does not depend on an a priori jet model. 
No "trigger particle" (parton proxy) is required and no 
"associated-particle" p t cuts are imposed. In the absence 
of trigger-associated pt cuts all minijet hadrons, which 
strongly overlap on p t those hadrons produced by soft 
processes (e.g., participant nucleon fragmentation along 
the collision axis), are accepted in the analysis. 

The phrase "minijet contribution" refers in the present 
context to the distribution of correlated hadron frag- 
ments from a minimum-bias parton energy spectrum av- 
eraged over a given A-A (or N-N) event ensemble. Be- 
cause the parton spectrum is rapidly falling (~ l/p?), 
with an observed lower bound near 3 GeV, the appar- 
ent minimum-bias parton spectrum is nearly monoencr- 
getic The term "minijets" then corresponds exper- 
imentally to jets localized near the 3 GeV lower bound 
(equivalent to parton en ergy scale Q ~ 6 GeV) , consistent 
with the original usage [26l \2l\. Minijets (minimum- bias 
jets) are further discussed in App. [S] 

In this analysis we report experimental tests of the 
local-thermalization hypothesis and conjectured bulk 
medium properties using minijets as probes of the sys- 
tem. By analogy with Brownian motion [28[ minijet 
probes (small-Q gluons) are just "large" enough (suf- 
ficiently energetic) to manifest as hadronic correlations 
(minijets) yet "small" enough to provide good sensitiv- 
ity to local medium properties and dynamics (e.g., other 
semihard partons) |29j |. 

It is essential to establish a theoretical baseline predic- 
tion for minijet correlations. In the absence of medium 
effects such correlations should correspond to a linear 
superposition of N-N collisions (binary collision scaling) 
as described by the Glauber model of A-A collisions 
(Glauber linear superposition or GLS). Minijets may be 
strongly modified in more-central collisions or even van- 
ish in an opaque thermalized medium [27, 30-32]. The 
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goal of this analysis is to determine where measured 
minijet correlations agree with baseline predictions (no 60 
medium effects) obtained from perturbative QCD as rep- 
resented by the hijing Monte Carlo [33| and to quantify 
any deviations from that baseline as a function of collision 
energy and centrality. Our results are further discussed 
in terms of the expected centrality trends for semihard 65 
parton scattering and fragmentation in dense, strongly- 
interacting media predicted by theoretical calculations 
and phenomcnological models. 

Angular correlations among the products from nuclear 
collisions are revealed by two-dimensional (2D) angular 
autocorrelations (Sec. [TTJ) defined on pseudorapidity and 
azimuth difference variables ?7a = ?7i ~ V2 and 4>a = 
4>i — 4>2 [H|-[36j]. Correlation sources include hadronic 
resonances, elliptic flow, quantum statistics (HBT) and 
semihard parton scattering (minijets). In proton-proton 
(p-p) collisions the observed an gula r correlations, when 
viewed using pair- wise pt cuts [37l l38j , are comprised 
of simple geometric structures: (i) a same-side (0a < 
7r/2) 2D peak at the origin on (^Aj^a), (ii) an away- 
side ridge in the form of dipole cos(0a — ?>")> and (iii) 
a ID peak on 77A centered at the origin, (i) and (ii), go 
with hadron p t > 0.35 GeV/c (for p-p collisions), are 
interpreted together as minijet angular correlations, and 
(mWalls mainly below hadron p t = 0.5 GeV/c [2(| H3, 

Other correlation analyses have been performed with 
RHIC data, but most have focused on specific features of 
angular correlations. Several PHENIX studies (e.g., (39|) 85 
were restricted to ID azimuth correlations. Other STAR 
and PHOBOS analyses have imposed so-called trigger- 
associated pt cuts (e.g., [i(|) which retain only part of 
the jet structure and reduce or exclude other contribu- 
tions. One other analysis [Il[ does consider p t -integral 90 
2D angular correlations (albeit over a restricted central- 
ity range) and is discussed further in Sec. IVIII Al 

The STAR Collaboration previously reported mea- 
surements of minimum-bias 2D angular correlations 
for charged-particle pairs from Au-Au collisions at 95 
130 GeV [2l|. Significant correlation structures from sev- 
eral sources were reported, including those interpreted as 
minijet contributions. Centrality variation of the same- 
side 2D peak was inconsistent with expectations from 
jet-quenching theory [27|, l3"oTj32l |. Instead of diminish- 100 
ing with increasing Au-Au centrality (as expected in jet 
quenching scenarios), the same-side peak amplitude in- 
creased strongly with centrality, and the azimuth width 
decreased instead of increasing. Most surprisingly, the 
width on relative pseudorapidity ?7a increased more than 
2-fold from peripheral to central collisions. However, thc 105 
limited statistics of the 130 GeV Au-Au data did not 
permit detailed study of the centrality dependence of the 
correlation structure. 

In the present analysis the method of Ref. [2l| has been 
applied to charged hadron production from minimum-no 
bias Au-Au collisions at -Vsnn = 62 and 200 GeV O . 



A preliminary report of results was presented in 43j . The 



much larger data volume (compared to the 130 GeV data) 
and two collision energies make possible a detailed study 
of the centrality and energy dependence of correlation 
systematics. The new results confirm our previous ob- 
servation of unexpected centrality trends [21], which in 
retrospect constitute the discovery of r\ broadening of the 
same-side peak, but also reveal for the first time the on- 
set of strong deviations from binary-collision scaling at a 
specific Au-Au centrality common to both energies. 

Taken together, our analysis results reveal that the cor- 
relation structure of interest (minijet structure) evolves 
with centrality according to a simple Glauber linear- 
superposition baseline, consistent with no novelty in A- 
A collisions compared to p-p, up to a specific centrality 
point where evolution of several parameters undergoes a 
sharp transition (large slope changes within a small cen- 
trality interval) to a qualitatively different smooth trend. 
The large increase in jet-like structure above the transi- 
tion point relative to the GLS trend contrasts with ex- 
pectations of strong jet quenching in more-central A-A 
collisions (26l . [29i - [3l| ] . The anomalous centrality evolu- 
tion then consists of the sharp transition and the un- 
expected increase in jet-like correlations in more-central 
Au-Au collisions relative to theoretical expectations, as 
discussed in Sees. IVII and IVIII TH 

Given the discovery of anomalous centrality evolu- 
tion involving correlations interpreted in p-p collisions in 
terms of minijets we wish to test various theoretical col- 
lision scenarios, especially those assuming rapid thermal- 
ization to form a dense bulk medium nominally opaque 
to jets. We hypothesize that pQCD minijet structure 
should follow a GLS reference trend in A-A collisions un- 
less modified by interactions leading to thcrmalization. 
We determine to what extent jet structure is modified 
from p-p to central Au-Au collisions relative to the theo- 
retical baseline. Perturbative QCD-based Monte Carlo 
model hijing [33| (without jet quenching) provides a 
nominal GLS theory baseline. In thermalization scenar- 
ios we expect to see strong reduction and other modifica- 
tions (symmetric peak broadening) or even extinction of 
jet-like correlations. If those expectations are not met we 
may question the theoretical assumptions. Recent corre 
lation predictions from the transport model AMPT 44 1 
and event-wise hydrodynamic model NEXSPHERIO 
are discussed in that context in Sec. IVIII .11 

This paper is organized as follows: The analysis 
method, data selection and measured angular correla- 
tions are described in Sees. HIllIVI respectively. The 
model-fitting procedure, fit results, and systematic un- 
certainties are presented in Sees. IWVTIl respectively. Re- 
sults and implications for heavy ion collision interpreta- 
tions are discussed in Sec. lVIlIl and a summary and con- 
clusions are presented in Sec. lIXI Further analysis details 
are presented in five Appendices. 
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II. ANALYSIS METHOD 

Number correlations (reported here) on binned two- 
particle momentum space (as opposed to p t correla- 
tions [34l . |36j]) are commonly reported as a ratio in each 
5 bin of the number of sibling pairs (from the same col- 
lision) to a number of reference or mixed pairs (from ! 
different but similar collisions - see App. [C]). The ratio 
relative to unity is 



Ap 

Piei 



Psib — Pmi 
Pref 



I: 



(1) 



where bin indices are suppressed, p S ib denotes the den- 
sity of sibling pairs, p rc f is the reference density of mixed 
pairs, and ratio r = p s ib/Vref- Expressions using binned 
pair counts explicitly are provided in App. [Cj The per- 
pair measure defined in Eq. ([T]) is useful for quantum 65 
correlations (46[, where e.g. the number of correlated 
pairs in a bin on invariant relative momentum may be 
approximately proportional to the single-particle density 
squared. However, correlation structures associated with 
initial-state scattering (relative to number of participant 70 
nucleons) or hadronization (relative to final-state hadron 
multiplicity) are better described by a per-particle mea- 
sure. Defined symbol Ap/^/p rc f represents such a per- 
particle measure designed specifically to test the null hy- 
pothesis that a nucleus-nucleus collision is equivalent to 75 
a Glauber linear superposition of N-N collisions. The 
statistical measure Ap/,/p rc f defined by 



Ap 



c ' t ^ = ^fZf(r-l) (2) 

Prcf v 



is equivalent to Pearson's normalized covariance (or cor- 
relation coefficient). The numerator Ap is the covariance 
of fluctuating particle numbers in two single-particle his- 
togram bins, and the denominator (effectively J p' rc{ ) is 
approximately the geometric mean of two single-particle 
number variances, leading (in the Poisson limit) to per- 
particle normalization. The explicit form is given in 
Eq. ([3]) and App.[Cj The right-most expression in Eq. ([2]) 
insures the cancelation of acceptance effects (and other 
experimental artifacts) in the ratio r. 

Prefactor -J p' le{ in Eq. @ is the ideal geometric-mean 
single-particle density absent acceptance and inefficiency 
effects, approximately the single-charged-particle density 
d 2 N c \ L /dr]d(j) averaged over the angular acceptance. Both 
62 and 200 GeV Au-Au multiplicity distributions are con- 
stant on pseudorapidity to within 1-2% for \-q\ < 1 (47i — 
|49| . The prefactor can therefore be approximated by 
y/ p rcf ~ N c h/ ArjAcf), where iV c h is the corrected mean' 
charged-particle multiplicity within the acceptance for 
each centrality bin (see Tables IIIII and IIV|) , and the an- 
gular acceptance for this analysis is defined by A77 = 2 
(\rj\ < 1) and A(j> = 2ir. The correlation measure used in 
this analysis is then 11 
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(Tab- I), (3) 



where CI denotes the charge-independent summation 
over four charge-pair combinations. r a f, is the sib- 
ling/mixed ratio of normalized total pair numbers in 2D 
histogram bin (a, b) averaged over charge-pair combina- 
tions, event-multiplicity bins (within a given centrality 
bin), and collision- vertex-position bins, the average being 
weighted by sibling-pair number as described in App. [Cl 

Indices (a, b) represent an unspecified 2D binning of 6D 
two-particle momentum space (pi,P2)- For p-p collisions 
it was shown that projections onto subspaces (pti,pt2) 
and (rjx, rj 2 , <j>\, <p2) are complementary (correlation struc- 
ture is factorized with minimal information loss) (37l . |38| . 
In this analysis, as in Ref. [2l| . only projections onto the 
angular subspace are reported (p t -integral correlations). 

In Ref. [35| correlation structures on 2D angular sub- 
spaces (rji, 772) and (<pi, 4>2) were found to be invariant on 
sum axes 771 +172 and 4>i+4>2 within the STAR TPC accep- 
tance. 4D angular subspace (771, 772, 4>i, $2) can then be 
simplified by projection onto difference axes 771 — 772 and 
4>i — 4>2 by averaging without loss of information along 
the sum axes within the TPC angular acceptance, thus 
forming a 2D angular autocorrelation [34l |36| |. Indices 
(a,b) then label 2D bins on difference axes (?7a,<M). 

The autocorrelation technique in the context of nuclear 
collisions applies to angular correlations only, and only 
in the case that invariance on the sum axes is a good ap- 
proximation (e.g., within restricted intervals on 77). The 
technique does not apply to correlations on (m t ,m t ) or 
(yt,yt) for instance. We use the formal term "autocorre- 
lation" initially for the purpose of definition and adopt 
the shorter form "correlations" subsequently in the text. 

This analysis is unique in part because it introduces 
several new techniques, including (a) consideration of the 
full range of A-A centralities down to N-N collisions, (b) 
application of a statistically well-defined per-particle cor- 
relation measure, (c) definition of a Glauber linear super- 
position reference, (d) accurate model fits to 2D angular 
correlations, (e) proper control of several systematic bi- 
ases including pile-up effects, distortions due to canonical 
suppression arising from centrality definition and distor- 
tions due to variation in position of the collision vertex 
and collision multiplicity. Those aspects are discussed 
further in Sec. IVIII A1 



III. DATA 

Data for this analysis were obtained with the STAR 
detector [42j using a 0.5 T uniform magnetic field paral- 
lel to the beam axis. Minimum-bias triggered events for 
collision energies v / snn = 62 and 200 GeV were obtained 
by requiring a coincidence of two Zero-Degree Calorime- 
ters (ZDCs) and a minimum number of charged-particle 
hits in the Central Trigger (scintillator) Barrel (CTB). 
Charged-particle measurements with the Time Projec- 
tion Chamber (TPC) and event triggering are described 
in |42j |. Primary vertices, defined using TPC tracks, were 
required to fall within 25 cm of the axial (z-axis) center of 



FIG. 1: (Color online) Perspective views of two-dimensional charge-independent angular correlations Ap/^/Pref on (7/a,</>a) 
for Au-Au collisions at ^/snn = 200 and 62 GeV (upper and lower rows respectively). Centrality increases left-to-right from 
most-peripheral to most-central. Corrected total cross-section fractions are (left to right) 84-93%, 55-64%, 18-28% and 0-5% 
for the 200 GeV data and 84-95%, 56-65%, 18-28% and 0-5% for the 62 GeV data (see Tables HTT1 and HV|) . 



the TPC. The data accepted for this analysis included 6.7 
million events at ^/snn = 62 GeV (Run 4 - 2004) and 1.2 35 
million events at 200 GeV (Run 2 - 2001). The present 
analysis is not limited by statistics; the 1.2M events from 
Run 2 are sufficient for all analysis requirements. 

Accepted particle trajectories fell within the TPC ac- 
ceptance defined by p t > 0.15 GeV/c, \rj\ < 1.0 and 2n m 
azimuth. Primary tracks in each event were required to 
have a Distance of Closest Approach (DCA) less than 
3 cm from the reconstructed primary vertex, accepting a 
large fraction of true primary hadrons plu s approximately 
12% background contamination [HI, |49| from weak de- 45 
cays and interactions with detector material. Conversion 
electron-positron backgrounds were reduced by excluding 
particles with dE/dx (specific energy loss in the TPC) 
within 1.5cr of that expected for electrons in the momen- 
tum ranges 0.2 < p < 0.45 GeV/c and 0.7 < p < 0.8 50 
GeV/c. Charge signs were determined, but particle iden- 
tification was not otherwise implemented. Further details 
of track definitions, efficiencies and quality cuts are de- 
scribed in Refs. [48l fHofl. 

Event pileup results in tracks from an untriggered 55 
event coexisting with a triggered event in the TPC. Al- 
though the pileup rate for these data was typically less 
than 1% such pileup can produce significant unwanted 
structure in angular correlations. A method to correct 
angular correlations for pileup is described in App. [D] 60 

The pileup-corrected minimum-bias event sample at 
each energy was divided into eleven centrality bins: nine 
each with nominally 10% of the total cross section and 
the most-central 10% split into 5% bins. The corrected 
centrality fractions reported in Tables lTHI and llVl were de- es 
termined from the minimum-bias distribution plotted as 



dN cvcnt /dN^ 4 versus N 1 ^ 4 on accepted event multiplic- 
ity N c h after adjustments for inefficiencies due to trigger- 
ing, collision vertex finding and particle trajectory recon- 
struction. That distribution is nearly uniform because 
the minimum-bias distribution dN event /dN c h is observed 
to approximate a "power-law" trend oc N~ } ^ 4 12]. The 

low-multiplicity end point of the distribution on N 1 ^ 4 
was constrained by measured p-p minimum-bias collision 
multiplicities [5l| normalized to the STAR TPC accep- 
tance. 

Multiplicity iV c h used to determine the centrality was 
integrated over the same pseudorapidity acceptance \r)\ < 
1 used for the correlation analysis. Use of iV c h from a re- 
stricted interval (e.g. \r)\ < 0.5 as in 18]) to define the 
collision centrality results in artifacts in 2D histograms 
due to canonical suppression. Correlations (fluctuations) 
within the restricted pair acceptance are suppressed rela- 
tive to those outside it, leading to substantial systematic 
errors in the angular correlations. 

Centrality is represented in a Glauber context by pa- 
rameter v = 2(iVbin) / (^Vpart)j the average number of N-N 
binary collisions per incident participant nucleon (in ei- 
ther nucleus) as obtained from Monte Carlo Glauber- 
model simulations [52J related to 62 and 200 GeV 
minimum-bias distributions on dN c h/dr] from [48], HU 
and denoted and ^200 • Parameter v is matched to 
observable N c h through the integrated total cross sec- 
tion via the approximately rectangular power-law dis- 
tribution on N 1 ^ 4 as described in Ref. [H} . At the 
lower-multiplicity end point (half-maximum point) v = 1 
while at the upper-multiplicity end point (6 = 0) v = 
5.29 ± 0.20 and 6.17 ± 0.23 for 62 and 200 GeV data, re- 
spectively. The estimated mean value of v for p-p (N-N) 



collisions is 1.25 (differing from 1 because of the skewness 
of the N-N multiplicity distribution [12|). 

The Glauber parameters can also be viewed as purely 
geometric measures unrelated to a particular N-N pro- 
cess: v can be thought of as the average participant path eo 
length. The 200 GeV parameters (assuming a 42 mb N-N 
cross section [Ejjj]) are then adopted as default geometry 
measures for both energies. Centrality measure v facili- 
tates tests of the N-N linear-superposition hypothesis. 

Estimates of ensemble-mean N c h for each centrality 65 
bin were obtained from minimum-bias multiplicity distri- 
butions [49| and from Monte Carlo Glauber-model sim- 
ulations assuming a two-component hadron production 
model The two methods agreed within 10% (most 
peripheral) and 1% (most central) and were within 6% ™ 
and 3% for the intervening centralities for the 200 and 
62 GeV data, respectively. Average values were used for 
the corrected multiplicities N c ^, listed as angular density 
dN ch /dr} = N ch /2 in Tables QID and El 



IV. 2D ANGULAR AUTOCORRELATIONS 

Figure Q] shows perspective views of data histograms 
Ap/ytprcfiVA, 4>a) for representative centralities ob- 80 
taincd from Au-Au collisions at v /snn = 62 and 200 GeV. 
The histograms show (within a constant factor) the 
event-wise mean number of correlated pairs per final- 
state particle in each (77^, 4>a) bin. 

The pair angular acceptances were divided into 25 85 
bins on the r]A axis and 25 bins on <j)\ 1 a compromise 
between statistical error magnitude and angular reso- 
lution. The histograms are by construction symmetric 
about ?7a = and 0a = 0, w. The 25 bins on </>a actu- 
ally span 27r + 7r/12 to insure centering of major peaks 
on azimuth bin centers. Statistical errors are ~ ±0.002 
(±0.004) for 62 (200) GeV data near \t)a\ = 0. Because of 
the rjA dependence of the p air acceptance s tatistical er- 90 
rors increase with \t)a\ as \J Arj/(Ar] — \t)a\) f° r V accep- 
tance A77 = 2 but are uniform on 4>a ■ Statistical errors 
are approximately independent of centrality for this per- 
particle measure. Statistical errors are larger than the 
above trends by approximately v2 for angle bins with 
?/A = 0, 4>a = and ±7r because of reflection symme- 
tries. An additional overall y/2 increase applies to the 
two most-central centrality bins which split the top 10% 95 
of the total cross section. 

Although the principal features of the correlations pre- 
sented and discussed in the remainder of this article are 
consistent with those reported in Ref. [2l| the details 
are much clearer. The centrality dependence is accu-100 
rately determined over the full range from N-N to b = 
Au-Au, and the collision-energy dependence is measured 
for the first time. The most-peripheral Au-Au central- 
ity bin corresponds approximately to minimum-bias N-N 
(~p-p) collisions. The null hypothesis that A-A colli-105 
sions are Glauber linear superpositions of N-N collisions 
is clearly valid for the more-peripheral Au-Au collisions, 



but strongly falsified for more-central collisions. 

Principal histogram features include (in the same order 
as Fig. [2] panels after the fit residuals): (i) a same-side 
(\4>a\ < 7r /2) 2D peak (approximately Gaussian) cen- 
tered at (t?Aj0a) = (0,0) which increases in amplitude, 
narrows on 4>a and dramatically broadens on t]a with 
increasing centrality; (ii) an ?7A-mdependent away-side 
(|0a| > t/2) dipole (ridge) clearly apparent for the most 
peripheral bin and strongly increasing in amplitude with 
increasing centrality; (iii) an ^-independent cos(2 0a) 
azimuth quadrupole with maximum amplitude for mid- 
central collisions. The quadrupole feature has been con- 
ventionally identified with elliptic flow, (iv) an approx- 
imately 0A-mdependent ID peak on rjA (approximately 
Gaussian) centered at ija — (observed along the front 
edge of some panels and/or superposed on the away-side 
ridge), diminishing in amplitude to zero with increasing 
centrality; and (v) a narrow 2D peak (approximately ex- 
ponential) at (t)a, 4>a) — (0, 0) (due mainly to conversion 
electrons and quantum correlations or HBT). 

The features observed in peripheral 200 GeV Au-Au 
collisions agree well with those reported previously for 
200 GeV p-p collisions [37],|38j]. Based on systematic stud- 
ies of two-particle angular and transverse-momentum 
correlations for p-p collisions [2(J H3, !H, [55[ we conclude 
that the same-side 2D peak [excluding the sharp spike 
at (0,0)] and away-side ridge represent semihard parton 
scattering and fragmentation (minijets). The visual fea- 
tures and fitting model components are discussed in the 
following section. 

V. MODEL FUNCTION AND 2D FITS TO DATA 

2D histograms have significant advantages over ID pro- 
jections and nongraphical numerical methods (e.g. some 
V2 analysis). Multiparameter fits to 2D histograms are 
generally less ambiguous than fits to their ID projections 
because covariances among fit parameters are reduced by 
the additional information in the 2D histograms. 



A. 2D model function 

2D angular correlation histograms from Au-Au colli- 
sions for 22 energy and centrality combinations were fit- 
ted with a six-component model function. The Au-Au 
model was adopted from one developed during analysis 
of 200 GeV p-p collisions [13, HI]. The fit model for p-p 
collisions was motivated by the simple geometrical forms 
apparent in the correlation data, not by an a priori phys- 
ical model. A cos(2 </>a) azimuth quadrupole component 
was added to the p-p model to describe the Au-Au data. 

The model function applied to Au-Au correlation his- 
tograms includes (in the same order as panels in Fig. [2] 
after the fit residuals) (a) a same-side (SS) 2D Gaussian 
on (va,4>a), (b) an 7/A-independent away-side (AS) az- 
imuth dipole cos(0a — tt), (c) an ^-independent azimuth 
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FIG. 2: (Color online) Fit decomposition of the 46-56% centrality data for 62 GeV Au-Au collisions. The upper panels show 
from left to right the corrected data, model fit, fit residuals (data — model) and same-side 2D Gaussian. The lower panels 
similarly show the away-side azimuth dipole, the nonjet azimuth quadrupole, the ID tja Gaussian and the 2D exponential. 
This centrality is just below the sharp transition at vtrans = 3. Fit residuals (c) are scaled up eight-fold relative to the data. 



quadrupole cos(2 0a), (d) a 0A-independent ID Gaus- 
sian on ?7a, (e) a narrow SS 2D exponential on {rj/\, 0a) 
and (f) a constant offset. The combined six-component 25 
model function in that order is 



F = At exp I -- 



<T<t>£ 



>1A 



' '/A 



(4) 



+ A D cos(0a - 71") 
+ Aq cos(2 a ) + A exp 

+ A 2 exp 




Given that mathematical description of the data the 
model elements can be interpreted physically. Terms 
(a) and (b) taken together are interpreted as a mini- 40 
jet contribution based on arguments in App. at least 
in peripheral Au-Au collisions. Term (c) is convention- 
ally identified with elliptic flow (5(| [57j- Term (d) is 
associated with participant-nucleon fragmentation (local 
charge conservation results in unlike-sign charged hadron 45 
pairs appearing nearby on r\ [37l HH). Term (e) mod- 
els quantum correlations (HBT) and conversion-electron 
pairs. 

The away-side ridge, attributed to pt conservation (e.g. 
back-to-back jets), can be modeled either by an AS az- 50 
imuth dipole (better for low-p t fragments from minimum- 
bias or small-Q partons) or by a ID Gaussian at 0a = tt 
with image peak at 0a = — tt (better for higher-p t frag- 
ments from more-energetic partons). With decreasing 
parton energy and increasing peak width the AS Gaus- 55 



sian periodic array approaches an AS azimuth dipole as 
a limiting case [59| . The AS dipole then provides a more 
efficient description of the AS ridge. The effect of the AS 
ridge model choice on other fit parameters is included in 
the systematic uncertainties discussed in Sec. VII. 

X 2 fits to the data were conducted by averaging the 
combined model function over a 5 x 5 grid within each 
(f?A, 0a) bin rather than using function values at bin mid- 
points. The averaging technique becomes important in 
regions where the model function has large curvatures. 
In particular it affects the relation between the 2D expo- 
nential and 2D Gaussian near the angular origin. 

Figure [2] shows an example of fit decomposition and 
residuals using the 62 GeV 46-56% corrected central- 
ity bin (nominal 50-60% bin). Similar results are ob- 
tained for each centrality bin and energy. The upper 
panels show (left to right) data, model fit, residuals 
(data— model) and SS 2D Gaussian. The lower pan- 
els show the AS azimuth dipole cos(0a — ""), azimuth 
quadrupole cos(20a), ID Gaussian on 77A and 2D expo- 
nential. For this centrality, and for all other data ex- 
cept a few more-central bins, the residuals are compara- 
ble in magnitude to statistical errors and are negligible 
compared to the amplitudes of the principal correlation 
structures. In this example the ID Gaussian on t?a (g) 
describes a small artifact in 2D correlations at tja = 

Absence of physically-significant structure in the fit 
residuals indicates that the 2D fit model of Eq. Q ex- 
hausts all statistical information in these data. The data 
do not require additional model components. Fit resid- 
uals for a few more-central bins at both energies include 
a small-amplitude non-statistical structure (AS dipole 
modulation on ?7a) discussed in Sec. lVII Ci| For minimum- 
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bias (pt-integral) angular correlations the SS 2D peak 
is well-described by a single 2D Gaussian. There is no 
systematically significant evidence for a separate non- 
Gaussian "ridge" in the SS 2D peak structure for angular 
5 correlations integrated over p t > O.f 5 GeV/c. Discussion eo 
of possible additional data structure and model compo- 
nents (e.g., V3) is presented in Sec. IVIlfHl and App. [Bl 



B. Model-fit results « 

Best-fit descriptions of data were based on a % 2 min- 

10 imization procedure. For most centralities any substan- 
tial excess contribution to total x 2 was confined to the 
acceptance edge \tja\ > 1.5. Excluding those bins from 70 
the fitting procedure had a negligible effect on the best-fit 
model parameter values. The resulting model parameters 

15 are presented in App . IFl ( Tables [TTT1 and HVf . The columns 
of Tables IIIII and IIVI correspond to the eleven centrality 
classes. The first eleven rows in both tables present the 
fit parameters from Eq. ^ plus the statistical (fitting) 
and systematic uncertainties. The remaining rows re- 

20 port centrality and other derived parameters. Centrality 75 
is measured by participant path length v from a 200 GeV 
Au-Au Monte Carlo Glauber model used as a common 
geometry parameter for both energies. Most of the model 
fit parameters exhibit strong variations with centrality. 

25 The error matrix for the fit parameters revealed statis- 80 
tically significant covariances among some of the param- 
eters, for example among the dipole, quadrupole and SS 
2D Gaussian amplitudes for the more-central histograms. 
In order to account for covariances the corresponding sta- 

30 tistical uncertainties were estimated by an iterative pro- 85 
cedure. A given parameter was displaced from its opti- 
mum x 2 fit value, the other ten parameters were adjusted 
to minimize x 2 > the selected parameter was further dis- 
placed and the data refit until the total x 2 increased by 

35 f. The reported uncertainties thus reflect covariances 90 
among the parameters. The incremental uncertainties 
(r.m.s. variances and covariances) are all observed to be 
small compared to the magnitudes of the parameters. 
Fitting errors for the model parameters in Eq. ^ are 

40 listed in App. 10 95 
The fit model in Eq. (H]) includes non-orthogonal com- 
ponents which could lead to ambiguities in the best-fit 
solutions (e.g., multiple local x 2 minima). The possibil- 
ity of ambiguities was studied in detail. Ambiguities were 

45 eliminated by conducting many independent x 2 fits as- 
suming thousands of initial- value combinations for the 11 
model parameters to locate the best-fit global minimum. 

For the more-central data at both energies a continu-100 
ous fitting ambiguity developed when the ID r\ Gaussian 

50 amplitude was allowed to become negative. The concave- 
upward shape in the away-side ridge for the more-central 
data (see Fig.Q]- right-most panels) pulled the ID Gaus- 
sian amplitude negative and forced the width to becomeios 
large. In combination with the ^-broadened SS 2D 

55 Gaussian that lead to a continuous fitting instability: 



The offset, dipole, quadrupole, 2D Gaussian amplitude 
and ID Gaussian amplitude and width could simultane- 
ously co-vary over substantial intervals without signifi- 
cantly reducing the residuals (always less than 5% of the 
SS 2D peak amplitude). 

The source of the ambiguity was identified as a 
statistically-significant residual structure not described 
by Eq. Qj. The ambiguity could be prevented by plac- 
ing a (lower or upper) bound on the value of any one of 
the affected model components. Since the m = 1, 2 sinu- 
soids and SS 2D peak parameters are the main focus of 
this paper we chose to remove the instability by requiring 
the ID r\ Gaussian to be non-negative. The impact of the 
imposed lower bound on the best-fit parameters was in- 
cluded in estimation of the systematic uncertainties. See 
Sec. IVII Cl for further discussion of the excess residuals. 



VI. ANOMALOUS CENTRALITY EVOLUTION 

Figure [3] shows the centrality and energy dependence of 
physically-relevant fit parameters reported in Tables IIIII 
and IIVI Two important trends emerge: (i) strong central- 
ity variation tightly correlated between the SS 2D peak 
and AS dipole amplitudes and (ii) smooth variation of 
the azimuth quadrupole amplitude. In this analysis we 
hypothesize that trend (i) is related to semihard parton 
scattering (minijets for more-peripheral collisions). The 
azimuth quadrupole [trend (ii)] is conventionally associ- 
ated with elliptic flow. Comparisons with previous Au- 
Au 130 GeV results are discussed in App. |EJ 

The term "anomalous" in the section title refers to two 
aspects of centrality evolution: (a) the large increase in 
slope of centrality trends for the SS 2D peak amplitude, 
its width on rj/\ (represented by model parameters A\ and 
a VA ) and AS ID peak amplitude Ad by factors 3.5, 5 and 
3.5 respectively within one centrality bin (sharp transi- 
tion) and (b) the large amplitude increase (up to twice 
the GLS trend) and significant azimuth width decrease 
of the SS 2D peak with increasing centrality, both trends 
contradicting conventional expectations for jet quenching 
in a strongly-coupled dense medium. Anomalous central- 
ity evolution is discussed further in Sec. IVIII Bl 



A. Centrality and energy trends 

With increasing centrality the SS 2D peak exhibits (a) 
a pronounced increase in the slope of the amplitude trend 
(i.e., a transition in the parameter trend with centrality) 
at transition point vtr&ns = 3.1 ± 0.3 (including statis- 
tical errors and bin-to-bin correlated and uncorrelated 
systematic uncertainties) accurately mirrored by the am- 
plitude trend of the AS dipole, (b) a similar increase in 
the slope of the 77A width at the same transition points 
vtrans-, and (c) a 4>a width decrease. There is no signif- 
icant difference in v tra ns for the two collision energies. 
<70 A for more-central collisions approaches a fixed value 
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FIG. 3: Fit parameters for (tjaj^a) correlation data from Au-Au collisions at ^/snn = 62 (open symbols) and 200 GeV 
(solid symbols) versus centrality measure v computed at fixed energy (200 GeV). The same-side 2D Gaussian amplitudes, t]a 
widths, and (j>A widths are shown in the left, center and right panels respectively of the upper row. The lower row shows from 
left to right the amplitudes for the dipole, quadrupole, and same-side peak width aspect ratio cr, ;A /cr^ A . Fitting errors are 
indicated by error bars where larger than the symbols. Solid lines connect the points for clarity. The dotted and dashed curves 
indicate Glauber linear superposition estimates for 62 and 200 GeV peak amplitudes respectively, as discussed in the text. The 
quadrupole data are consistent with Ref. [6(j. The hatched regions indicate the full range of systematic uncertainties listed in 
App. [F] The vertical dark bands indicate estimated v equivalents for N-N collisions and b = Au-Au collisions. 



fa 0.7. Above the transition point the SS 2D peak and 
AS dipole amplitude trends for both energies increase 
uniformly on centrality to v pa 5, beyond which they de- 
crease. The correlated-pair yield decreases above v = 5 25 
are intriguing, but are also comparable to the systematic 
uncertainties presented in App. |FJ 

The SS 2D peak is actually strongly elongated on az- 
imuth (o0 a :c, ;a = 2:1) in peripheral collisions. But with 
increasing centrality the angular asymmetry reverses and 30 
the 2D Gaussian becomes three times broader on t]a than 
on </>a- The smooth shape evolution is shown by the as- 
pect ratio plotted in Fig. |3] (bottom-right panel). It is 
notable that the SS 2D peaks for 55-64% 200 GeV and 
56-65% 62 GeV histograms in Fig.[T]have unit aspect ra- 35 
tio (equal r.m.s. widths on jja and <?!>a), but the peaks 
appear to be elongated on 77A because the histograms as 
plotted have an aspect ratio of 27r:4 ps 3:2. The SS 2D 
peak widths are further discussed in Sec. I VIII CI 

In contrast to the sharp transition in same-side 2D 40 
peak properties, the azimuth quadrupole amplitude 
varies smoothly with centrality, with no manifestation 



of the transition behavior observed in the SS 2D peak 
trends. The quadrupole amplitude depends only on geo- 
metric path length v (estimated by 1^200) , with functional 
form independent of collision energy [60] . 

The energy dependence of the SS 2D peak in num- 
ber angular correlations can be compared with that of 
the azimuth quadrupole and the SS peak in previously- 
measured pt angular correlations [3g|. In Ref. [6(| an 
inferred energy factor of the form ln(y / SNN _ /13.5 GeV) 
was found to describe t>2 data measured by ^4q = 
2poO) v!{2D}(&) (denning 2D fit parameter v 2 {2D}) 
above 17 GeV, where po(b) — dN c h/2Trdr] is the single- 
particle 2D angular density. The quadrupole amplitudes 
obtained in this analysis agree with those from Ref. [6(| ■ 

In the present analysis we observe that the same- 
side 2D peak amplitudes for two energies and central 
Au-Au collisions are in the ratio A\ (62) jA\ (200) = 0.57± 
0.06 (stat.) which can be compared with the energy- 
factor ratio ln(62.4/13.5 GeV)/ ln(200/13.5 GeV) = 0.57. 
The energy dependence of the SS peak amplitude in num- 
ber correlations is also consistent (within systematic un- 
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certainties) with a \n(y/stm/10 GeV) energy dependence 
of the SS peak amplitude inferred from p t angular corre- 
lations in Ref. |36| . Given the uncertainties in the lower- 
energy SPS pt correlation measurements 10 GeV can be 
interpreted as a lower limit on the intercept consistent 55 
with 13.5 GeV from Ref. [g^. Thus, the SS 2D peak and 
azimuth quadrupole collision-energy trends agree above 
17 GeV and depend only on 1h(v/5nn)' The latter de- 
pendence is consistent with QCD processes. 

The ID peak on ?ja> interpreted to arise from eo 
participant- nucleon fragmentation [H, [6l| , is small com- 
pared to the SS 2D peak and falls monotonically to zero 
by mid-centrality [y ~ 3.5). 

Centrality parameter v = 2Nbin/N par t is smoothly 
(not discontinuously) related to the fractional cross sec- 65 
tion cr{b)/(T tot and to participant number -/Vp art . The 
sharp transition in SS 2D peak properties near v = 3 
appears as well when fit parameters arc plotted on other 
centrality measures. In the case of V part the transition 
shifts to the extreme left end of the parameter range and 70 
is therefore visually obscured. Parameter v presents the 
essential linear-superposition reference in a simple form: 
proportionality of the reference to binary-collision num- 
ber iVbin as discussed in the next subsection. 



B. Testing the linear-superposition hypothesis 

Accurate measurement of centrality trends for Au-Au 
angular correlations down to the N-N limiting case 1 
makes possible a rigorous comparison of Au-Au correla- 
tions to N-N binary-collision scaling — the Glauber linear- 
superposition reference (baseline). In the GLS refer- 
ence model of Au-Au collisions the SS 2D peak ampli- 
tudes (and volumes) from minimum-bias p-p (~N-N) col- 1 
lisions are linearly superposed (summed) at the angular- 
difference origin (77 a = 4>a = 0) proportional to the 
Glauber- model number of N-N binary collisions N [ n . In 
the GLS hypothesis the SS 2D peak 77A and 0a widths 
retain fixed values characteristic of p-p collisions. 

For per-particle measure A/j/^/p rc f in Eq. ([3]) binary- 
collision scaling of the SS 2D peak amplitude and vol- 
ume translates to scaling as ^Vbi n /-V c h- If X pp represents , 
a correlation peak amplitude or volume in p-p (~N-N) 
collisions the GLS variation with A-A centrality should 
be 



Xaa(v) = X p 



X 



Abin A c h, p 

A c h,AA 
V 



vN, 



ch, pp 



(2/A part )Af ch ,AA 



i'i> 



l + x(v- 1)' 



(5) 



where the second line assumes the two-component hadron 
production model of Kharzeev and Nardi (K-N) [54j . 
Amplitude or volume X pp can be estimated by direct p-p 
measurements or by extrapolation to N-N from several 100 
peripheral A-A centralities. 

Parameter x, the coefficient of the binary-collision scal- 
ing component, is held fixed in the K-N two-component 



model. Assuming x to be independent of centrality pro- 
vides a reasonable description of experimental probabil- 
ity distributions on multiplicity 52]. More-differential 
spectrum analysis suggests that the effective x increases 
substantially from p-p to central Au-Au collisions (62i — 
lo4j |. For this GLS reference x is held fixed at 0.02, the 
p-p value for acceptance A77 = 2 [25|, [63| . 

GLS references for the SS 2D and AS peak amplitudes 
are shown as the dotted and dashed curves in Fig. [3] for 
62 and 200 GeV data respectively. The amplitude data 
closely follow the GLS reference with increasing central- 
ity (within small systematic uncertainties) until the tran- 
sition point I'txansi beyond which the data substantially 
exceed the reference trends. Peak widths on both 77 and 
<t> show significant deviations from GLS constant values 
°7ja = 0-55 and =1.10 corresponding to N-N (p-p) 
collisions. SS peak width trends are further discussed in 
Sec. IVIII Gl 

The aspect ratio trend in the lower-right panel is par- 
ticularly interesting. It confirms the large eccentricity 
of the same-side 2D peak observed previously in p-p 
collisions with substantial elongation on <p (<T0 A :<t ??a = 
2:1) [38| and shows the strong evolution with Au-Au cen- 
trality to large elongation on 77 (o' VA :a ( j >A = 3:1). In p-p 
collisions the elongation on (f> was found to vary strongly 
with particle p t , with larger <f) elongation for smaller par- 
ticle pt down to 0.5 GeV/c for each particle [381 ] - 

HIJING [33| predictions for the SS 2D peak ampli- 
tude [43| from 200 GeV Au-Au collisions with jets imple- 
mented but no jet quenching deviate strongly in more- 
central Au-Au collisions from the GLS trend extrapo- 
lated from p-p data. A discussion of the discrepancy is 
presented in Sec. IVIII II The hijing SS 2D peak widths 
on (va,4>a) are respectively 0.75 and 0.9 (radians) and 
remain constant with centrality, in marked contrast to 
the large angular asymmetries and strong centrality de- 
pendence observed in the data. 



VII. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties in the parameters of the fit- 
ting function in Eq. ((4]) are primarily due to secondary 
particle and other contamination backgrounds in the 
data, uncorrected detector and event reconstruction ef- 
fects, ambiguities in the choice of fitting model function, 
and statistically significant residual structures not ac- 
counted for by the fitting model. The specific sources 
of uncertainty and the method of error estimation are 
discussed in the following subsections. Systematic uncer- 
tainties for the fitting parameters are listed in Tables IIIII 
and HV| in App. El 



A. Uncertainties in the histogrammed data 

The largest source of systematic uncertainty is a 12% 
non-primary particle contamination (48l . |49| with un- 
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known correlation structure in the particle sample used 
for the analysis. This background is primarily from 60 
weak-decay daughters from the collision and secondary 
particles produced in the detector material which were 
misidentified as primary particles, i.e. those emitted di- 
rectly from the triggered collision. e + -e~ pair contam- 
ination produced by photon conversions in the detector 65 
material are discussed in the following subsection. Cor- 
relation measure Ap/ y/p re { was computed assuming par- 
ticle DCA < 3 cm (standard cut admitting a 12% sec- 
ondary contamination) and DCA < 1 cm (reduced con- 
tamination fraction) and the resulting histog rams were 70 
compared. Any difference in correlation structure should 
be dominated by secondary particles preferentially re- 
moved by the modified DCA cut. Differences were found 
to be dominated by statistical fluctuations. Any system- 
atic structures were less than 3% of the primary corre- 75 
lation amplitudes, resulting in a ±3% uncertainty esti- 
mate assigned to the five amplitude parameters in the fit 
model. 

Pileup contamination was corrected as described in 
App. |D] We observe that pileup mainly affects the ID so 
77 Gaussian amplitude and mainly near mid-centrality 
91], causing the amplitude to vary non-monotonically 
with centrality v. Comparing the centrality dependences 
of parameter Aq before and after pileup correction sug- 
gests that about ±15% of the full pileup effect may 
remain in the Au-Au 62 GeV data (with larger ini- 
tial pileup fraction) after correction. Residual pileup 
contamination in the 62 and 200 GeV correlations was 

85 

therefore estimated separately in each centrality bin as 
±0.15Ap/7p^(pileup) from Eq. (IDT) . 

Pair reconstruction inefficiencies [21( induce depletion 
of Apj y'/Oref at small opening angles, visible in more- 
central collisions as grooves in uncorrected Ap/ y/p rc f 90 
near (t7A)</>a) = (0,0) for \t)a\ < 0.08 and |0 A | < 1. 
Although corrections (pair cuts to both sibling and ref- 
erence pairs) remove most of this effect, close examina- 
tion of the 2D histograms suggests that small artifacts 
remain which are approximated as a 2D Gaussian with 95 

amplitude 0.025 and -0.04 for 62 and 200 GeV 0-5% 

centrality data respectively, and with 77A and 0a widths 
0.08 and 0.5. Estimates for the other centrality bins were 
obtained by scaling the above amplitudes by the particle 
pair density, (d 2 N ch / drjdfi) 2 . 100 

Other systematic effects considered include: intermit- 
tent electronics outages, pseudorapidity acceptance de- 
pendence on longitudinal (z-axis) collision vertex posi- 
tion in the TPC, collision- vertex position inaccuracy due 
to reconstruction error, particle momentum resolution, 105 
TPC central-membrane particle-trajectory crossing inef- 
ficiency, and residual dependence on the event-mixing bin 
sizes for collision vertex position in the TPC and event 
multiplicity. The overall contribution to 2D angular cor- 
relations from those sources was found to be insignificant 110 
compared to the reported correlation structure. 

The effect of the above uncertainties in 2D correlation 
histograms on the fit-model parameters was estimated by 



separately adding each of the above representations [e.g. 
0. 15 Ap/y/p IC { (pileup), and small Gaussians for errors due 
to two-particle inefficiency, electronics outages and pseu- 
dorapidity acceptance dependence] to the data, refitting 
the data, increasing the amplitude of the added function, 
refitting the data again, and so on until a linear trend ex- 
ceeding statistical fluctuations could be determined. 

Parameter uncertainties due to secondary backgrounds 
were assumed to be Gaussian distributed. Those due 
to residual pileup were assumed to be uniformly dis- 
tributed, implying that the "true" parameter value lies 
between fitted values obtained by adding or subtract- 
ing 0. 15 Ap/y/p^f (pileup) to the 2D correlations with 
uniform probability. The other parameter uncertainties 
were assumed to be uni-directional with uniform proba- 
bility Mean shifts and variances estimated in this 
subsection were added linearly to those discussed in the 
next subsection. 

The uncertainty in pre-factor J p' Ie{ or d 2 N c \ 1 /dr]d(j) de- 
rived from spectrum analysis is ±8% (±7%) for the 62 
(200) GeV data OBI. This overall normalization un- 
certainty is not included in the parameter uncertainties 
reported in Tables Hill and HVl in App. [Fj 



B. Uncertainties arising from the fitting model 

In Sec. [V] we noted that the away-side ridge can be 
described by either an azimuth dipole or periodic ar- 
ray of ID Gaussians with common width ct^as (59j - In 
addition, for more-central collisions the away-side ridge 
amplitude is not constant on 77^ but displays a concave- 
upward dependence, a feature not readily described by 
the fitting function adopted for this analysis. Also, simi- 
lar model fits applied to 2D angular correlation data with 
pt cuts imposed suggest that the SS 2D peak for those 
data may be better described by a non-Gaussian func- 
tion. Alternative model descriptions of the SS 2D peak 
for more-central collisions could include a SS ID Gaus- 
sian ridge on azimuth (see Sec. IVIimi and App. IBl. We 
therefore estimate the extent to which the choice of fitting 
model function affects the accuracy of our description of 
the principal correlation structures. 

To explore systematic uncertainties derived from the 
choice of model function the components in Eq. (U) were 
modified. The data were refit and any changes in the pa- 
rameters of the unmodified components in Eq. (j4]) were 
recorded. The modifications included: a periodic series 
of ID Gaussians replacing the AS azimuth dipole, addi- 
tional i] 2 ^ cos (f>A and 77^ cos 2(/a terms (modeling alter- 
native ?7a dependence), modified SS 2D Gaussian [dif- 
ference exponent n allowed to deviate from 2 in e.g. 
exp{— (x — a;)"}], description of the SS 2D peak for more- 
central collisions as the sum of a ID Gaussian on <fi& 
plus an alternative 2D Gaussian, modified 2D exponen- 
tial (difference exponent allowed to deviate from 1), and 
similar exponent variation for the ID 77 Gaussian. Cor- 
responding shifts in the parameters of unaltered compo- 
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FIG. 4: (Color online) (a) 2D data histogram from 0-5% central 200 GeV Au-Au collisions [(0, 0) bin suppressed], (b) Additional 
model component fitted to the 2D histogram and consisting of an AS azimuth dipole modulated by function F(t]a) described in 
the text, (c) Fit residuals including the new model component, (d) Data histogram in the first panel minus the additional model 
component and fitted offset. Remaining structure is described accurately by an away-side dipole and same-side 2D Gaussian. 



nents of Eq. @ for each modified fit model determined 40 
the uncertainties, assumed to be unidirectional with uni- 
form probability. Variances were obtained from the pa- 
rameter shifts. 

The sharp spike at (0,0) apparent in Fig. [5] is predom- 
inantly caused by e + -e~ pairs produced by photon in- 45 
teractions in detector material which survive particle- 
identification (dE/dx) and primary-particle selection 
cuts. Quantum correlations also contribute near the an- 
gular origin, as shown by projecting HBT correlations 
onto (f?A,<M)- Both background correlations are well 50 
described by the single 2D exponential in Eq. In 
addition to allowing the difference exponent to vary we 
also studied the impact of this component by remov- 
ing the 2D exponential from the model and a few his- 
togram bins near (0,0) from the fit and refitting the re- 55 
maining data with a truncated 8-parameter model. The 
resulting changes in parameters were assumed to be uni- 
directional, with uniform probability distributions. 



C. Additional model element for central collisions 

2D histograms for the most-central Au-Au collisions 
exhibit a distinct 77A dependence in the AS dipole not 
observed in less-central Au-Au or p-p collisions, with a 65 
minimum at t]a — 0. As an example, the data histogram 
for 0-5% central 200 GeV Au-Au collisions is shown in 
Fig. S] (first panel). Attempts to model the visible AS 
minimum with a Gaussian negative on 7/a and uniform 
on </>a were rejected because the fits were unstable and 70 
did not significantly reduce the residuals. 

Detailed examination of the 2D residuals (data — 
model) for the nominal fits in Sec.[V]for the most-central 
collisions at both energies led to modification of the AS 
dipole model component shown in Fig. 0] (second panel) : 75 
adding an AS dipole times an f^A-symmetric function 
with minimum value zero at the origin having the form 
F(va) = \VA\ m exp{-|W^72}, with m ~ 2, n ~ 5 
and a ~ 1.5. The fractional ryA modulation of the total 
AS dipole is about 15% for 0-5% central collisions and de- ao 
creases to zero for mid-central collisions. The additional 



model component is obviously orthogonal to the azimuth 
quadrupole component and any other v m for m > 2. As 
with other 2D model elements the AS dipole modulation 
was introduced in response to observed data structure 
and is not motivated by a physical mechanism. 

Residuals from a x 2 fit with the additional model com- 
ponent are shown in Fig. 0] (third panel) . Subtracting the 
AS dipole modulation term from the data histogram in 
the first panel leads to the fourth panel. What remains 
is a uniform AS dipole and an SS 2D peak well described 
by a 2D Gaussian. Those — 5% data are consistent with 
zero quadrupole amplitude and with the corresponding 
entries in Table IIIII This result suggests that in some 
cases apparent deviations of the SS 2D peak from an ideal 
2D Gaussian shape may actually result from superposi- 
tion of a small AS dipole modulation, with characteristic 
azimuth width (~ 1.6) much larger than the azimuth 
width of the SS 2D peak (cr 0A « 0.65). The SS 2D peak 
itself then does not deviate significantly from a 2D Gaus- 
sian. 

The impact of this residual on parametrizations of 
the most central collisions at 62 and 200 GeV was es- 
timated in two ways. First, the ID Gaussian compo- 
nent in Eq. (U) was replaced by the above residual model 
function F(j)a) cos (0a — tt), the data were re- fitted, and 
the shifts in the parameters of the remaining components 
recorded. Errors for the other centralities were estimated 
by scaling the preceding 0—5% errors by the approximate 
overall amplitudes of the residual in the centrality bin. 
Residuals above the transition centrality have shapes 
similar to the above model while those below the tran- 
sition are dominated by statistical fluctuations. Second, 
the residual histograms at each centrality were separately 
fitted with a dipole, quadrupole, SS 2D Gaussian and 
ID Gaussian on t]a- The resulting very small sinusoid 
amplitudes were included in the errors. The associated 
Gaussians were added to the nominal fitted Gaussians, 
and effective amplitudes and widths were computed from 
the volumes and second moments of the combined Gaus- 
sians. The differences between the effective amplitudes 
and widths and the nominal parameters were included 
in the systematic uncertainties. Each of the systematic 
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uncertainties in this subsection was assumed to be uni- 
directional, with uniform probability distribution. 



D. Total systematic uncertainties 

The mean shifts (in Sees. IVII Al and IVII CI) and total 
variances of the systematic uncertainty contributions (15 
in all) discussed in this section were summed for each of 
the eleven parameters in the model function, accounting 
for the symmetric or uni-directional nature of the uncer- 
tainty and its assumed probability distribution, Gaussian 
or uniform. If the mean shift exceeded the total sys- 
tematic uncertainty r.m.s. value, then the error bar was 
extended to include the nominal fit value. 

The principle sources of systematic uncertainty var- 
ied with centrality and parameter, but in general the 
secondary-particle contamination, the apparent r]A de- 
pendence of the AS dipole and the amplitude and shape 
(exponent) of the 2D exponential dominated. System- 
atic uncertainties due to some sources of systematic error 
tend to be correlated across parts of the centrality range. 
Total systematic uncertainties in adjacent centrality bins 
may therefore be partially correlated. The nominal fit 
parameter values and their statistical and systematic un- 
certainties are listed in Tables ITIT1 and ITVl in App. |FJ 

The total systematic uncertainties are also represented 
in Fig. [3] by the hatched regions at the bottoms of the 
panels. The full range of uncertainty is represented. In 
most cases the uncertainties tend to be symmetric about 
the plotted values. The exception is 200 GeV cr„ A where 
the uncertainties extend mainly above the plotted values. 



AS dipole should thus be relatively insignificant in more- 
central collisions. Additional energy- and momentum- 
conservation-induced correlations (e.g. p z \ x p z2 [66|]) 
would produce distinct tja dependence (hyperbolic func- 
tions) in the 2D angular correlations which are evidently 
too small to be statistically significant in the residuals. 

Any reduction of v 2 with increasing \rj\ should pro- 
duce a corresponding reduction in the quadrupole am- 
plitude with | 77a |- The systematic uncertainty estimated 
via the rj A cos 2<M model component allows for this pos- 
sible structure in the correlations. The V2{rj) data [67j 
for 1 77 1 < 1 do not require such a reduction, but would be 
consistent with reductions of a few percent. The present 
data do not require an | 77 a I -dependent quadrupole ampli- 
tude. Directed flow (vi) [68|, [69| might contribute an |tia| 
dependence to the dipole amplitude, but is estimated to 
be too small to observe in these data and is not apparent 
in the data or residuals. Higher-order azimuth sinusoid 
components were found to be negligible in the fit residu- 
als projected onto 4>A compared to the additional dipole 
component discussed in Sec. IVII CI Also, see Sec. IVIIlTIl 

It has been conjectured that a sextupole component 
referred to as V3 may actually be present in angular cor- 
relation data. Alternatively, the SS 2D peak may include 
a SS ID Gaussian component uniform on tia- Such struc- 
tures are not observed in fit residuals from p t -intcgral 
data using the standard model defined by Eq. The 
consequences of adding such model elements to the stan- 
dard fit model are discussed in App. [Bj In App. IB 51 re- 
lated systematic changes in the fit results reported in Ta- 
bles QTI] and QV] are estimated to be substantially smaller 
than the combined errors reported in those tables. 



E. Other possible correlation structures 

85 

Small variations near |?7a| = 2 in the cos(?i</>a) sinu- 
soids result from residual effects of finite collision- vertex- 
position and event-multiplicity bin size. Random struc- 
tures near \rj&\ = 2 reflect limited two- particle statistics 
near the 77 acceptance boundary. These structures are 9° 
found to have a negligible effect on fitted parameters. 

Simulated angular correlation structures due to reso- 
nance decays (mainly p° and lu for the present pt ac- 
ceptance) were found to contribute less than 10% of 
the 2D same-side Gaussian peak within |?7a| < 0.5 and 95 
I'M I < 2 [U [65| (and were negligible elsewhere). Such 
correlation structures were not observed in the fit resid- 
uals and thus were not included in the fitting model. 

Global transverse-momentum conservation produces 
per-pair angular correlations measured by Ap/p re f pro- 
portional to pn -pt2/N ch = {pnPt2/N ch ) cos(</>a) [Sa| and 
is therefore included in the AS dipole amplitude. The 
magnitude of the corresponding per-particle dipole am-100 
plitude in Ap/y/p lc f could be as large as 0.015 to 0.02 
at either energy, but should be independent of centrality. 
Any global momentum-conservation contribution to the 



VIII. DISCUSSION 

The sharp transition in jet-like angular correlation 
trends [93j revealed by the present analysis introduces 
a surprising new aspect of RHIC Au-Au collisions. The 
transition occurs at a specific value of mean partici- 
pant path length v common to both collision energies 
(y trans ~ 3). Such a transition in the SS 2D peak ampli- 
tude and width on r\A could mark the onset of a new cor- 
relation mechanism beyond semihard parton scattering 
and fragmentation. However, any proposed theoretical 
description, including novel collision mechanisms, must 
describe accurately the smooth centrality dependence of 
the SS 2D peak azimuth width and pt angular correla- 
tions [59] above and below t'trans- We now consider fur- 
ther details. 



A. Goals and unique aspects of this analysis 

This analysis has as its primary goal accurate descrip- 
tion of 2D angular correlations over the complete range 
of Au-Au centralities with a minimal complement of sim- 
ple functional forms. We wish to determine under what 
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FIG. 5: Fit parameters for (tja, 4>a) correlation data from Au-Au collisions at ^/snn = 62 (open symbols) and 200 GeV (solid 
symbols) versus centrality measure v illustrating a sharp transition in centrality trends at Vtrans ~ 3. The dashed curves labeled 
GLS are calculated predictions corresponding to linear superposition of N-N collisions (binary collision scaling). The bold solid 
lines illustrate slope changes by factors 3.5 in panels (a) and (b) and 5 in panel (c) within one centrality bin on v. The dotted 
line in panel (b) simply continues the slope trend from below the transition. Panel (a) shows that 62 and 200 GeV same-side 
2D peak amplitudes A\ are related by a log(y / ijvjv) factor ~ 0.6 [721 ]. whereas panel (c) demonstrates that the away-side ID 
peak amplitude Ad is approximately independent of collision energy. 



circumstances and to what extent more-central Au-Au 
collisions deviate from a simple linear superposition of 
N-N collisions according to the Glauber model of A-A 
collisions extrapolated from peripheral collisions. A ma- 
jor issue for this study is the extent to which minimum- « 
bias jet structure observed in p-p collisions is modified 
in more-central Au-Au collisions. Based on related stud- 
ies we adopt a minijet hypothesis: that minimum-bias 
jet-like structure in more-peripheral A-A collisions corre- 
sponds approximately to that from 3 GeV jets (minijets). 45 
We then examine the extent to which that structure is 
modified in more-central Au-Au collisions. The extent 
and nature of such modifications is compared to several 
theoretical scenarios as a test of their validity. 

The present analysis is unique in several aspects: (a) 50 
consideration of the full range of A-A centralities down to 
N-N collisions, (b) use of a statistically well-defined per- 
particle correlation measure, (c) definition of a Glauber 
linear superposition reference, (d) accurate model fits to 
2D angular correlations and (e) distinct measurements 55 
of a same-side 2D peak, away-side ID peak and nonjet 
quadrupole. Conventional V2 analysis [70] describes only 
ID azimuth projections and only with one Fourier se- 
ries term, disregarding the r\ dependence critical for dis- 
tinguishing separate structures and mechanisms. Analy- 
sis of "dihadron" azimuth correlations also considers 
only ID azimuth projections and subtracts a background 
based on possibly-biased V2 data f5il . Trigger-associated 
ID and 2D angular correlations [7l| include only part of 
the jet structure (high-pt hadrons) and typically do not 
represent the 2D structure with model functions. 

Those characteristics can be compared with a 2D an- 
gular correlation analysis by the PHOBOS collaboration 65 
presented in Ref. 4l|. In the PHOBOS analysis corre- 
lations were measured over 27r azimuth and r\ € [—3,3]. 



The correlation measure defined in Eq. (1) of Ref. [4J| 
is that used in Ref. [21]. Particle pt could not be recon- 
structed over the full (??, (^acceptance and was not used 
for the analysis in Ref. [41]. Although 2D angular cor- 
relations were inferred, emphasis was placed on the pro- 
jected ID distribution on ?7a which has a single peak fea- 
ture ("short range" correlations) fitted by a ID Gaussian 
(cluster model) . No distinction was made between the SS 
2D peak and ID Gaussian on t?a isolated in the present 
analysis. The underlying physical mechanism for the pro- 
jected ID peak was not specified. Properties of the fitted 
ID Gaussian were inferred only for the upper 50% of the 
total cross section due to systematic uncertainties in the 
event selection procedure for low-multiplicity events. The 
sharp transition in SS 2D peak properties reported in the 
present study is then inaccessible. It was concluded that 
angular correlation structure in central Au-Au collisions 
is similar to that in p-p collisions, and an inferred clus- 
ter size decreases with increasing centrality, dramatically 
contradicting the large increases in jet-like structure ob- 
served in the present study. 



B. Anomalous evolution of correlation structure 

The reported anomalous evolution has two aspects (see 
Fig. 3): (a) Three correlation model parameters (SS peak 
amplitude, AS peak amplitude, SS peak rj width) un- 
dergo large slope changes in their centrality trends within 
a small interval on centrality at v = ftrans common to 
two energies. In addition, the azimuth width of the SS 
peak, which decreases significantly from p-p up to ^trans, 
maintains a fixed value above that point, (b) The large 
increase in the SS peak amplitude above t'trans, inter- 
preted as a minijet manifestation, is inconsistent with 
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expectations of scattered-parton thermalization (strong 
jet quenching) in more-central A-A collisions, and the w 
decrease in the azimuth width is inconsistent with par- 
ton multiple scattering and other jet-quenching scenar- 
ios [13, Sa HI, H ■ 

Figure [5] illustrates manifestation of anomalous cen- 
trality evolution in the form of feature (a) large slope 
change within a small centrality interval. In all three 
panels the sharp transition (rapid slope change) occurs 65 
near Vt ians = 3 (slope change of solid lines). In the first 
two panels the slope ratios of the solid lines are fixed 
at 3.5. In the third panel the slope ratio is 5. Data 
deviations from the solid lines in the first two panels are 
consistent with fit uncertainties (error bars) . In the third 
panel the 62 GeV transition may be slightly displaced to m 
larger v. Nevertheless, three fit parameters at two en- 
ergies (six instances) exhibit large slope changes in their 
centrality trends within the same small centrality inter- 
val ftrans = 3.1 ± 0.3. The interval 0.6 about v = 3 
corresponds to 10% on fractional cross section a/ao- n 

In the first panel the 62 GeV data are brought into 
coincidence with the 200 GeV data by a factor 1/0.63. 
That relation corresponds to the energy scaling quantity 
R{y/sWN = 62 GeV) = log(62/13.5)/log(200/13.5) = 
0.57 from Ref. [72] which describes the azimuth so 
quadrupole per-particle energy dependence (relative to 
values at 200 GeV) reported there. The correspondence 
is notable. Also, whereas the SS 2D peak amplitude does 
scale as log(y / s]v r /v) the AS dipole amplitude in panel (c) 
does not. That difference is expected for pQCD dijets, ss 
holds for most Au-Au centralities, and the same energy 
scaling trend continues up to LHC energies [73l |. 

Centrality variation of the SS 2D peak amplitude and r\ 
width in panel (b) indicates that the integrated number of 
SS peak correlated pairs per final-state charged particle 90 
(SS peak volume) exceeds binary-collision scaling by an 
order of magnitude in more-central Au-Au collisions. In 
the minijet context the same-side peak pair number cor- 
responds to the product of the event- wise minijet number 
in the angular acceptance and the mean fragment pair 
number (« mean jet fragment multiplicity squared) [62j . 
Is the large correlated-pair increase due to excess produc- 
tion of minijcts with N-N properties in central A-A col- 
lisions (relative to binary-collision scaling)? Or does the 
mean multiplicity associated with each minijet increase™ 
relative to that for N-N collisions, the minijet number 
remaining consistent with binary-collision scaling? The 
quantitative correspondence among correlations, spectra 
and pQCD discussed in Sec. IVIIIEl seems to indicate that 
parton scattering changes little with increasing centrality, 
but the details of parton fragmentation to jets changes 
substantially. 

While the sharp transition in SS 2D peak properties isios 
itself notable, the fact that minijet correlations increase 
at all with increasing Au-Au centrality seems to con- 
flict with the conventional expectation that most jets are 
"quenched" in the dense medium formed in central A-A 
collisions, therefore not appearing as correlation struc-uo 



tures in the final state. Evidence from Raa measure- 
ments (hadron suppression at high p t ) [H| and high-p t 
jet AS azimuth correlations (disappearance of the away- 
side jet) seemed consistent with that expectation. In 
contrast, we observe that pt-integral minimum-bias jet- 
like structure increases dramatically with centrality. 



C. Anomalous SS 2D peak width trends 

The most notable feature of the SS 2D peak in more- 
central Au-Au collisions, its large elongation on ?7a, is 
not predicted by present pQCD theory. But neither is 
the comparable elongation on </>a observed in p-p colli- 
sions. A possible ?7a elongation mechanism arising from 
color connections between struck partons and their par- 
ent nucleons has recently been suggested [25|. The in- 
terplay between the SS peak widths is here described in 
more detail. 

Figure [5] shows centrality evolution of the two angular 
widths of the SS 2D peak in Au-Au collisions. The third 
panel shows the width (or aspect) ratios. The cr^ A trend 
(first panel) is approximated below the transition point 
at v « 3 by er, ;A = 0.53 + 0.2(j/ — 1) (dash-dotted line). 
The aspect ratio trends (third panel, dashed and dotted 
lines) are described over a larger centrality interval by 
(T ??A /o'0 A = exp{(v — vq)/1.8}, with vq = 2.25, and 2.55 
for 200 and 62 GeV respectively. The dashed and dotted 
curves describing (i0 A below the transition point in the 
middle panel are simply derived from those two results 
and indicate the consistency of the description. 

Nominal GLS trends assuming linear superposition 
of p-p collisions are indicated by the hatched regions. 
The predicted widths from HIJING are a VA = 0.75 and 
00 A = 0.9 (radians) independent of centrality, in marked 
contrast with the large angular asymmetries and strong 
centrality dependence of the SS 2D peak observed in the 
data. While the SS and AS peak amplitudes follow the 
GLS trend below the transition point (see Fig. 3) the SS 
peak widths do not. The individual width trends (the 
slopes) change substantially at the transition point, but 
the aspect ratio varies smoothly (exponentially) from p-p 
to more-central Au-Au collisions, exhibiting no sign of a 
slope change. Any viable theoretical description of angu- 
lar correlations in Au-Au collisions must accommodate 
that complex phenomenology. 



D. Minijets and pt angular correlations 

Are new (e.g., non-pQCD) collision mechanisms re- 
quired to accommodate the observed anomalous central- 
ity evolution? Information from p t angular correlations 
and single-particle spectra may help to reduce the am- 
biguity, pt angular correlations, complementary to num- 
ber angular correlations from the present analysis, have 
been obtained by inversion of ipt) fluctuation scale de- 
pendence from 200 GeV Au-Au collisions [H [M Izl- 
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FIG. 6: Same-side 2D peak width trends for (jja^a) correlation data from Au-Au collisions at ^/snn = 62 (open symbols) 
and 200 GeV (solid symbols) versus centrality measure v computed at fixed energy (200 GeV). The same-side 2D Gaussian j?a 
and (j>A widths are shown in panels (a) and (b) respectively. Panel (c) shows the same-side peak width aspect ratio o",, A /<70 A . 
The hatched regions indicate GLS trends based on p-p collision data. The bold curves are explained in the text. 



The pt correlation structure is qualitatively similar to 
that presented here (e.g., same-side 2D peak, away-side 
ridge, quadrupole), but there are significant quantita- 
tive differences. The same-side 2D peak amplitude trend 
on centrality, while increasing at least as fast as binary- 
collision scaling until v w 4 where it starts to decrease, 
does not show a substantial change in slope at f trans- 40 

In the minijct context p t angular correlations suggest 
that N-N semihard parton scattering continues to drive 
the SS 2D peak structure above u trans where the number 
of minijets (parton scatters) increases at least as fast as 
iVbin through and above i^rans until v 4. In this pic- 45 
ture it follows that the large increase of SS peak pairs 
in Fig. |3] is then due to strong modification of parton 
fragmentation leading to a large increase in the number 
of jet-correlated hadrons per initial-state parton scatter. 
That conclusion is substantiated by comparisons among 50 
spectrum hard components, pQCD fragment distribu- 
tions and minijet number angular correlations pBI . [62[ • 



E. Minijets and single-particle spectra 

The large changes in the number of correlated parti- 
cle pairs observed in this analysis should also be mani- 
fested in single-particle pt spectra. In Ref. 



a differ 



ential two-component analysis was applied to unidenti- 
fied charged hadron p t spectra from non-single-diffractive 
(NSD) p-p collisions at y/s = 200 GeV. Hard and soft 
spectrum components were isolated as limiting cases of so 
spectrum evolution with event multiplicity (centrality). 
The soft spectrum component was estimated from the 
iich — > spectrum limit. Subtracting the soft compo- 
nent from the full spectrum revealed the hard compo- 
nent centrality evolution. The hard component is inter- 65 
preted as hadron fragments from minimum-bias large- 
angle-scattered partons (minijets). The soft component 
is interpreted as fragments from dissociation of projectile 



nucleons. The hard component has been subsequently 
identified with jet-like p-p correlations [37l [3S| . 



In Ref. |6J] a similar two-component analysis was ap- 
plied to 200 GeV Au-Au spectra for identified pions and 
protons. Hard and soft spectrum components were iden- 
tified for each species. The hard components evolve 
strongly with centrality, accounting for essentially all of 
the per-participant spectrum evolution. Strong suppres- 
sion at larger p t is accompanied by much larger enhance- 
ment at smaller p t relative to binary-collision scaling of 
the N-N hard components. Suppression (at highp t ) and 
enhancement (at low pt) variations with centrality are 
strongly correlated, implying the same underlying mech- 
anism. The sharp transition on centrality for each hadron 
species matches the transition revealed in the present cor- 
relation analysis. 



In Ref. [25| the spectra from Ref. [64| were compared 
with a pQCD calculation of fragment distributions based 
on measured jet fragmentation functions (LEP, CDF) 
and pQCD predicted parton spectra. Again the agree- 
ment was found to be very good, lending strong support 
to interpretation of both spectrum hard components and 
minimum-bias jet-like correlations as pQCD jets. 

In Ref. [(32} preliminary data from the present analy- 
sis were combined with a pQCD prediction of jet num- 
ber in A-A collisions to infer parton fragment yields cor- 
responding to minijct production. The fragment yields 
were in turn compared with yields inferred from p t spec- 
tra for identified hadrons and the agreement was found 
to be good [si]]. From that exercise it was concluded 
that about one third of the final state in 200 GeV central 
Au-Au collisions is contained in resolved minijets, mainly 
from 3 GeV scattered partons. 
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F. Interpretation of the away-side dipole 

In Fig. [3] (left panels) the amplitudes of the away-side 
dipole structure and same-side 2D Gaussian follow the 
same centrality trend, strongly suggesting that they share 
a common mechanism. The mean energy ~ 3 GeV for 
minimum-bias scattered partons (2jJ_75| is comparable 
to the mean intrinsic kt ~ 1 GeV/cjTq within projectile 
nucleons, implying large acoplanarities for semihard scat- 
tered parton pairs and a broad away-side azimuth "ridge" 
(back-to-back parton correlation). It is easy to demon- 
strate that the large width of the away-side peak (~ 
Gaussian on azimuth at <fi — n) plus peak periodicit y o n 
azimuth are equivalent to dipole trend cos(</)a — tt) [591 ]. 

An away-side dipole could also be produced by global 
transverse-momentum conservation (Sec. IVIIEl and [66|) 
which might account for part of the observed AS dipole 
amplitude in peripheral Au-Au (N-N) collisions. For the 
per-pair quantity Ap/p re f that contribution to the dipole 
term would be proportional to l/N^, For per-particle 
measure A/?/^/p r cf used in this analysis that contribution 
should be independent of centrality. We observe that the 
AS dipole amplitude follows the same strongly-increasing 
centrality trend as the SS 2D peak. The close correspon- 
dence implies that the AS dipole is indeed a manifestation 
of transverse-momentum conservation, but at the parton- 
parton scale (dijets), not the nucleus-nucleus scale. 
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FIG. 7: Same-side 2D peaks from 38-46% central (left panel) 
and 9-18% central (right panel) 200 GeV Au-Au collisions. 
These histograms were obtained by subtracting the AS dipole 
(Ad), quadrupole (Aq) and offset (^3) elements of the 2D 
model fits from the data histograms. 




G. Other correlation structures 

The azimuth quadrupole component cos(2 0a) has 
been conventionally identified with elliptic flow. The 2D 
quadrupole amplitude reported here is related to con- 
ventional measure V2 by A Q = 2p a (b)v 2 1 {2D}, where 
Po(b) = dN c h/2irdr] is the single-particle angular density 
and "2D" denotes inference of quadrupole amplitudes 
from model fits to 2D angular correlations as described 
in this article [H, (77[. In Fig. [3] (bottom center) the 
smooth centrality variation of the quadrupole amplitude 
is in marked contrast to the SS peak amplitude and rj 55 
width. The quadrupole amplitude shows no counterpart 
to the sharp transition in SS peak properties. Although 
the centrality trends for SS peak properties and azimuth 
quadrupole are very different, the two amplitudes, when 
measured by statistically equivalent quantities Aq and «> 
Ai, share the same 1h.(^/snn) dependence characteristic 
of QCD scattering processes, as discussed in Sec. I VI AI 
Detailed analysis of the quadrupole component and its 
relation to other V2 methods is presented in Refs. [6011721]. 

The ID peak on 77A, interpreted to represent longitu- 
dinal projectile-nucleon fragmentation, has a simple cen- 
trality dependence. It is quite visible in p-p (N-N) colli- 
sions, especially for low-p t particles, |37j but falls mono- 65 
tonically to zero by mid-central A-A collisions (See Fig. [T] 
and Tables Mil and IIVI) . The structure was first observed 
at the ISR 61] and has been compared to similar struc- 
ture predicted by the Lund model (PYTHIA) [H,[78j|. In 



FIG. 8: Left: 2D histogram from Fig. [4] (fourth panel) pro- 
jected onto (/>a (light histogram). The bold histogram with 
near-zero amplitude is a projection of 2D residuals from Fig. [4] 
(third panel). The dash-dotted and dashed curves represent 
components of the model fit from the present analysis. Right: 
The bold histogram is the ID data histogram in the left panel 
minus the fitted AS dipole component (dashed curve). The 
dash-dotted curve is the ID projection of the fitted 2D Gaus- 
sian from this analysis. The bold dotted curves represent 
multipoles m = 1 ... 4 derived from the SS ID Gaussian (see 
the text). 



Ref. [78| a comparison of charge-dependent (CD) struc- 
ture dominated by that model component (Aq) was com- 
pared to PYTHIA data (Figs. 3 and 4 of that reference), 
and good agreement was observed. In contrast to Aq 
(nonjet quadrupole) and A\ (minijets), the Aq amplitude 
shows negligible energy dependence. 



H. Presence of higher multipoles measured by v m 

There are two modes in which higher-order (m > 2) az- 
imuth multipoles may arise in model descriptions of the 
present data. First, the data may require additional mul- 
tipoles for a statistically satisfactory description. That 
topic is mentioned below and discussed in more detail 
in App. |B] Second, the "standard" 2D model compo- 
nents defined in Eq. (j4]) may be decomposed into tja- 
independent multipoles (79j . based either on ID projec- 
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tions onto azimuth or on the full 2D model elements. 
The former is discussed in this subsection, the latter in 
App.E 

Figure [7] shows 2D angular correlations for two Au- 
Au centralities. The AS dipole (m = 1) and quadrupole 
(to = 2) terms of the 2D fit model have been subtracted 
from the data histograms. All that remains is the SS 2D 65 
peak, described by a 2D Gaussian. The 38-46% centrality 
(left panel) is just above the sharp transition at v trans = 
3. The 9-18% centrality (right panel) is at v = 5, the 
maximum of the SS 2D peak amplitude trend. 

As noted in Sec. IV Al the present fit model exhausts 70 
all statistical information in most combinations of cen- 
trality and energy. In more-central collisions significant 
residuals structure does remain, but only associated with 
the AS dipole (m = 1) as discussed in Sec. I VII CI Thus, 
any higher multipoles v m for m > 2 must come from 
the SS 2D peak structure shown in Fig. since the con- 
jectured v m are orthogonal to the subtracted model el- 
ements (to = 1, 2). However, the SS 2D peak has a 
strong curvature on r)/± which cannot be described by a 
ID Fourier series on azimuth. Thus, the SS 2D Gaussian 
in pi-integral angular correlations is unique. 

In this argument we focus on the 0-5% centrality bin 
where the Aq amplitude is consistent with zero, provid- 
ing a particularly simple example. However, the same 
arguments concerning higher v m apply to all Au-Au cen- ^ 
tralities, as demonstrated in Fig. [Jj 

Figure [5] (left panel) shows the 2D histogram data 
in Fig. [4] (fourth panel) projected onto azimuth differ- 
ence 4>a (light histogram). The dash-dotted and dashed 
curves represent the SS ID Gaussian (projected SS 2D 
Gaussian) and AS dipole respectively, derived from the 
fit to 2D histogram data as reported in Table Mil The 
bold histogram with near-zero amplitude is the ID pro- 
jection of the residuals in Fig. @] (third panel). The r.m.s. 
residuals amplitude is about 0.5% of the SS peak ampli- 
tude and consistent with statistical uncertainties. Thus, 95 
a model function consisting of SS ID Gaussian (two pa- 
rameters Aid, °>a) P ms AS dipole (one parameter An) 
exhausts all statistical information in the ID data his- 
togram. There is no necessity for additional Fourier com- 
ponents to represent these 2D angular correlations, as 100 
shown in Fig. [5] (left panel). That is, the systematic un- 
certainty interval for such amplitudes includes zero. 

Figure [5] (right panel) shows the data histogram in 
the left panel minus the fitted AS dipole (dashed curve). 
The difference (bold histogram) is described to the sta- 105 
tistical limits of data by a ID Gaussian (dash-dotted 
curve) with amplitude Aid, consistent with the resid- 
uals in the left panel. The bold dotted sinusoids (SS- 
peak multipoles) in the right panel represent the first 
four Fourier components (to = 1 ... 4) of the SS 1D 110 
Gaussian measured by 2po(b)v^ n = F^o^^) Am, with 
F m (cr<f, A ) = cr0 A \/277rexp(-TO 2 cr|^/2). The to = 2 
term estimates the "nonflow" contribution from the SS 
2D peak to the total azimuth quadrupole. The to = 3 
term is the azimuth sextupole ("triangular flow" [80j),ii5 



the to = 4 term (just visible) is the octupole component. 

For the 0-5% central data in Fig. [5] we obtain Am = 
0.585 ± 0.06, o> A = 0.63, F 3 (0.63) = 0.077 and 
Po(0-5%) = 107, from which we predict U3 = 0.015 ± 
0.0008. We also obtain v 2 = 0.025 ± 0.0013 (nonflow) 
and Vi — 0.007 ± 0.0004. The uncertainty estimates are 
based on the « 10% total uncertainty in A\ which is 
the dominant source. Similar results are obtained for 
other centralities, but the nonjet quadrupole measured 
by u|{2D} then contributes substantially. In general, at- 
tempts to measure a u 3 (or higher multipole) component 
in ID projections onto azimuth of 2D angular correla- 
tions can be anticipated accurately from the properties 
of the SS 2D peak determined by this analysis. 

The 2D structure of angular correlations imposes 
strong constraints on proposed data models. For the 
p t -integral data presented in this analysis the SS 2D 
Gaussian is the most efficient description of same-side an- 
gular correlations (requiring the fewest independent pa- 
rameters). Any Fourier-series representation of the SS 
2D peak (or hybrid representation including both 2D 
Gaussian and Fourier components) would be less effi- 
cient, requiring more parameters. Whereas the nonjet 
quadrupole measured by u|{2D} arj d obtained from fits 
to 2D histograms has negligible curvature (is approxi- 
mately uniform) on tja within the STAR TPC acceptance 
any multipoles associated with the SS 2D peak must all 
exhibit the same large curvature corresponding to the to\ 
dependence of the SS 2D peak. 

Since 2D angular correlations projected onto ID az- 
imuth implicitly include a W3 Fourier component, among 
others, as part of the SS 2D Gaussian, the SS 2D peak 
contribution to U3 and any other inferred Fourier com- 
ponent must be acknowledged before any claim of an in- 
dependent v m component is made. Any higher-u m con- 
jecture should be presented in the context of the full 2D 
histograms, not just ID projections onto azimuth. Inter- 
pretation of ID Fourier components of the SS 2D peak 
as flow manifestations competes with jet production by 
fragmentation as an alternative mechanism. From this 
analysis we conclude that a sextupole Fourier component 
t>3 is not required by the 2D data. 

Although the t] interval external to the STAR TPC 
acceptance is obviously excluded from this analysis, any 
correlation structure there should be considered in rela- 
tion to 2D model fits within the TPC acceptance as well 
as ID Fourier series. 

The observation of anomalous centrality trends re- 
ported in this paper applies to the SS 2D peak in the 
correlation data, as represented by the required model ele- 
ments defined in Eq. ^ . For alternative model represen- 
tations, for example those which invoke higher azimuth 
multipoles, the relevant parameter combinations remain 
those that describe the amplitude and widths of the SS 
2D peak in the data, regardless of how that peak struc- 
ture may be represented mathematically. See App. [Blfor 
further discussion. 



20 



< 0.8 



0.6 



7 (a) 



0.4 



0.2 




L(b) 



Au-Au data 




HIJING 



GLS ^ °- 25 



^NN.p-p 

'///, 



— 0.45 



0.4^ (c) 

'-■-?///?/////&//////,&■////////£. 
HIJING 



6 

V 



-| 0.2^ 
c a 0.15^ 

3 o.i 

05 




; ^ 0.05 W//. 



p-p data 




HIJING 



6 

V 



6 

V 



FIG. 9: (a) Amplitude of SS 2D peak (jet-correlated pairs per final-state hadron) for data (solid points), GLS extrapolation 
from p-p data (dashed curve) and from the hijing Monte Carlo (open points, dotted curve), (b) Single-particle production 
extrapolated from p-p data (dashed curve, GLS), and from more-central Au-Au data (solid line) and hijing (dash-dotted line, 
open points), (c) Single-particle hard-component production per N-N binary collision from p-p data, from more-central Au-Au 
data and from hijing (open points), (d) Jet-correlated pair production (SS 2D peak) per N-N binary collision extrapolated 
from p-p data (GLS), and from Au-Au data (solid points) and hijing (open points). 



I. Centrality trends for hijing 

40 

The hijing Monte Carlo with jet quenching off is nom- 
inally a linear superposition of N-N collisions modeled by 
PYTHIA within the context of a Glauber model of A-A col- 
lisions, hijing should therefore provide a GLS reference 
for Au-Au collisions but fails to do so. In this subsection 45 
we discuss the discrepancy. 

Figure [5] (first panel) shows SS 2D peak amplitude A\ 
for default hijing vl.382 [H[ simulations of 200 GeV 
Au-Au collisions (open points) compared to 200 GeV 
data (solid points) from Fig. [3] (first panel ), hijing is 50 
specifically formulated to describe minijets [3j]. If hard 
scattering is disabled in hijing the SS 2D peak disap- 
pears (Ai -> 0). The GLS (N-N) extrapolation from 
p-p data [dashed curve, see Eq. (JSJ)] is A\ — 0.065^/[l + 
0.02(t/ — 1)]. The dotted curve (explained in the text 55 
below) is 0.05f/[l + 0.16(f - 1)]. Apparently, hijing de- 
viates strongly from the GLS reference extrapolated from 
p-p data, whereas we expect hijing (no jet quenching) 
to represent an equivalent N-N linear superposition. The 
source of the discrepancy is determined as follows. eo 

Figure [5] (second panel) shows single-particle yields 
from hijing (points) and those expected from a GLS 
superposition of N-N in Au-Au collisions inferred from 
p-p spectra (63|. The p-p extrapolated (GLS) per- 
participant-pair multiplicity in K-N two-component 65 
form [5l] is 2.5[1 + 0.02(i/— 1)]. The hijing equivalent is 
2.1 [1 + 0.16(1/ — 1)]. We find a large difference in hard- 
component coefficient x: 0.02 vs 0.16. Also, hijing soft 
component Snn is 20% lower than p-p data. The K-N 
trend for more-central Au-Au data is 2.5[1 + 0.1(i/ — 1)]. 70 

Figure [5] (third panel) shows the hard-component (jet- 
related) hadron angular density per N-N binary collision 
based on p-p data, inferred from more-central Au-Au 
data and extracted from hijing (open points) in the 
second panel. Again we find that the hijing hard- 75 
component yield per N-N binary collision is much larger 
than what is observed in p-p collisions and even larger 



than observed in central Au-Au collisions. 

Part of the explanation for the hijing discrepancy 
comes from the underlying parton spectrum. The default 
parton spectrum cutoff (lower bound) in hijing is po = 2 
GeV/c. The spectrum cutoff inferred from 200 GeV p-p 
spectrum data is 3 GeV/c [25j, |63j. That difference has 
substantial consequences. The parton spectrum varies 
approximately as p^ 6 near 3 GeV/c. A shift in the lower 
bound by factor 2 /3 leads to an increased hard-scattering 
cross section per N-N collision by factor (3/2) 6 ps 10. 
We observe in Fig. [5] (third panel) that the hijing hard- 
component single-particle yield is a factor 7 larger than 
that inferred from p-p data, consistent with the difference 
in effective spectrum cutoffs. Another part of the hijing 
explanation comes from correlated-pair yields. 

Figure [5] (fourth panel) shows the 2D den- 
sity of jet-correlated pairs per N-N binary collision 
(l/Nbi n )p (b)Ai, with p (b)Ai = Ap[SS peak] expressed 
in terms of an absolute pair- number density and po(b) = 
dN c h/2irdri pa y/p Ic i, the single-particle density. What is 
plotted then is a jet-related pair-number density per N- 
N binary collision. We find that the jet-correlated pair 
density from hijing is about 2/3 that observed in p-p 
collisions, not seven times larger. From the pair densi- 
ties in this plot we calculate coefficients X pp for reference 
curves in the first panel: 0.05 = 2ir x 0.017/2.1 (hijing) 
and 0.065 = 2n x 0.026/2.5 (GLS trend from p-p data). 

We conclude that whereas hijing binary-collision- 
scaled single-particle yields are much greater than those 
for p-p (or even Au-Au) data, consistent with its default 
parton spectrum cutoff, the corresponding number of jet- 
correlated particle pairs is substantially lower than for 
p-p data. Spectrum hard components and jet angular 
correlations extracted from measured p-p data are mutu- 
ally consistent with measured parton fragmentation sys- 
tematics and a pQCD parton spectrum with cutoff near 
3 GeV [25|, [62|, [63[ . When combined they accurately pre- 
dict the centrality trend for more-peripheral Au-Au colli- 
sions (GLS extrapolation of p-p data), hijing spectrum 
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and correlations hard components are mutually consis- 
tent only if a very different parton fragmentation scheme 
is associated with the default 2 GeV parton spectrum 
cutoff. The combination does not predict the central- 
ity trend for more-peripheral Au-Au collisions (deviates w 
strongly from GLS extrapolation of p-p data) . 

Because of the parton spectrum power-law trend most 
scattered partons appear near the spectrum lower bound, 
which for default hijing is 2 GeV. The parton fragmen- 
tation process is poorly known there (see Ref. [25} for 65 
absence of accurate fragmentation function data below 
about 7 GeV parton energy). Fragmentation to charged- 
hadron pairs which could contribute to the SS jet peak, 
as modeled by hijing, appears to be substantially less 
than for real N-N collisions. The apparent agreement 70 
between hijing and minijet correlation data near v = 1 
seems to be an accident: hijing underestimates both p-p 
spectrum soft-component Snn and the number of jet- 
correlated pairs, the errors nearly canceling in the ratio. 
With increasing Au-Au centrality excess hard-component 75 
single-particle production in hijing combines with too 
few jet-correlated pairs to produce the observed large de- 
viations from the GLS trend extrapolated from p-p data. 
Implementing the hijing jet quenching mode causes re- 
duction of the same-side 2D peak amplitude with increas- ao 
ing centrality, further deviating from the Au-Au data. 



J. Centrality trends for other conventional models 

85 

We compare centrality trends of Au-Au angular corre- 
lations from the present analysis to those expected for 
three generic models of nucleus- nucleus collisions: a) 
Glauber linear superposition of N-N collisions, b) par- 
ton/hadron rescattering in a dissipative medium, and c) 90 
a locally-thermalized "opaque" medium. 

a. Linear superposition of N-N collisions — If 
Au-Au collisions were simply linear superpositions of 
N-N collisions (A-A transparency) then parton scattering 
and jet production in N-N collisions should be randomly 95 
superposed in A-A momentum space according to the 
Glauber-model number of binary collisions. The same- 
side 2D peak widths should remain constant with cen- 
trality while the amplitude should increase according to 
Eq. ([5]) . The away-side ridge should remain fixed in shape™ 
and follow the amplitude of the same-side peak. The 
shape of the single-particle spectrum hard component 
should also remain unchanged, and its amplitude should 
follow binary-collision scaling [64[. The GLS model is 
in fact our baseline reference for the centrality evolu-105 
tion of A-A angular correlations and single-particle spec- 
tra. It describes minijet correlations and spectra well 
for v < ^trans, but fails dramatically above the transi- 
tion and does not describe the azimuth quadrupolc. It is 
worth noting that the quadrupole amplitude reaches atuo 
least half its maximum value within the centrality range 
where N-N linear superposition accurately describes the 
SS 2D peak and dipole (Fig.©. 



b. Parton/hadron rescattering in a medium — Par- 
ton and/or hadron rescattering in a dissipative medium 
(e.g. as in cascade or transport models) implies that a pri- 
mary scattered parton and its hadron (jet) fragments are 
randomly deflected in angle and lose energy /momentum 
to the medium, as in Brownian motion [28| . The r\ and 
<f> widths of hadron-number and pt same-side 2D peaks 
should both increase. The same-side 2D peak amplitude 
for number correlations may be reduced compared to the 
binary-collision reference. The pt-correlations peak am- 
plitude should definitely decrease, and the spectrum hard 
component should be shifted to lower momentum and 
possibly reduced in amplitude (parton dissipation, "jet 
quenching" ) . The data reveal that the r\ width of the SS 
2D peak increases, but the azimuth width decreases with 
increasing centrality, inconsistent with random multiple 
scattering. A reduction in the p t same-side peak am- 
plitude is observed in data [74], but not until v > 4, 
well above t'trans- And the same-side peak amplitude for 
number-correlations increases dramatically above ftrans- 

Calculations using the multi-phase transport Monte 
Carlo model AMPT [44j for trigger-associated particle dis- 
tributions on 7?a and <pA from central Pb-Pb collisions at 
y/SNN — 2.76 TeV and central Au-Au collisions at y/s^N 
= 200 GeV were presented in Refs. [8l| and [82j respec- 
tively. Both papers include "string melting," or a par- 
ton cascade, as well as fluctuating initial conditions and 
pQCD jets described by hijing. With Lund-model [H| 
fragmentation parameters and parton-cascade cross sec- 
tions adjusted to fit V2 values at mid-centrality, the model 
predicts an 77-broadened SS 2D peak and a broad double- 
peaked AS ID structure on azimuth. The authors at- 
tribute those structures to anisotropic flow from fluctuat- 
ing initial conditions (initial-state densities) and medium 
response to jets but do not quantitatively distinguish the 
two sources. The predicted double-peaked AS structure 
is not observed in the present data, and neither is the pre- 
dicted double-hump (on n^) SS peak in Ref. 82j. AMPT 
predictions for the measurement conditions reported in 
this analysis are not available. 

c. Thermalized opaque medium — In some models 
of A-A collisions copious initial-state low-energy scat- 
tered partons are expected to contribute to formation 
of a locall y-th ermalized opaque medium with zero mean 
free path [2TL l3(il - [32l | . In those models low-energy scat- 
tered parton thermalization is assumed to occur almost 
immediately after impact (< 1 fm/c), and the scattered- 
parton energy is converted into thermal energy within 1- 
4 fm/c [3l], 32 1 which raises the local temperature ("hot 
spots"), while the parton momentum, being conserved 
overall, is dissipated to many final-state hadrons whose 
angular distributions arc broadened. With increasing 
centrality same-side number and pt angular correlations 
from low-energy partons (minijets) should be broadened 
and greatly reduced in amplitude relative to the GLS ref- 
erence, as estimated in [29j. The opposite trends are ob- 
served for the SS 2D peak amplitude and azimuth width 
while the r)A broadening reported here far exceeds that 
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expected from early-stage minijet interactions [29l . l3ll .l32T] 
and jet quenching (33[. 

Predictions for per-pair (Ap/p re f) angular correlations 
from Au-Au collisions at ^snn — 200 GeV for charged 
particles with 0.2 < p t < 2.0 GeV/c from the non- viscous 6 o 
3+1 hydrodynamics event- wise model NEXSPHERIO [45[ 
for four centralities from a 0-80% cross section fraction 
are presented in Ref. (83|. The nexspherio model as- 
sumes that initial conditions are described by the nexus 
event generator [84j which includes event- wise initial- 6 5 
state density fluctuations and soft, semi-hard and hard 
scattering processes. The generated 2D angular corre- 
lations reveal a SS 2D peak and AS ID ridge, but the 
physical origin of those structures is not identified [83| . 
Dramatic increases in the SS 2D peak amplitude and ??a ™ 
width in more-central Au-Au collisions reported in the 
present analysis are not predicted by the model. In- 
stead, nexspherio predicts that the r]A width should 
decrease. Differential study and identification of phys- 
ical mechanisms for correlation structures predicted by 75 
the nexspherio model are anticipated. 



IX. SUMMARY AND CONCLUSIONS 

80 

We have measured charged-particle number angular 
correlations from Au-Au collisions at >/snn = 62 and 
200 GeV projected onto relative azimuth </>a and relative 
pseudorapidity tja for eleven centrality bins on 0-95% of 
the total Au-Au cross section. The dominant features 85 
are a same-side 2D peak (approximately Gaussian), an 
away-side cos(^a — t) dipole and a cos(2</>a) quadrupole. 
The same-side 2D peak and away-side dipole in more- 
peripheral A-A collisions, first observed as correlation 
structures in minimum-bias p-p collisions and directly 90 
related to features of p-p single-particle spectra, can be 
reasonably interpreted as corresponding to minijets from 
minimum-bias parton scattering and fragmentation. 

p-p spectrum and correlation data combined with 
assumed Glauber-model linear superposition (binary- 95 
collision scaling) provide an essential baseline reference 
for heavy-ion collisions. The amplitude of the same-side 
2D peak in Au-Au collisions approximately follows N-N 
binary-collision scaling from peripheral to mid-central 
collisions and the widths are slowly varying, extending™ 
the p-p minijet hypothesis across the centrality range 
v < ftrans in heavy ion collisions. 

It is conventional practice in analysis of RHIC data 
to omit a substantial fraction of the complete A-A cen- 
trality range from observation (i.e., more-peripheral col-">5 
lisions). The essential GLS reference is then inaccessi- 
ble, and it is not clear what is unique about A-A col- 
lisions relative to p-p or N-N collisions, or where ex- 
ceptional behavior actually appears on centrality. The 
present analysis establishes for the first time that jet-no 
related correlation data from Au-Au collisions actually 
closely follow the GLS trend (binary-collision scaling) up 
to mid-centrality. Strong deviations of 2D angular cor- 



relations from the GLS reference in more-central colli- 
sions are quantitatively determined relative to the GLS 
reference. Such information provides a more differential 
(hence rigorous) test for any proposed theoretical model. 

At an intermediate centrality denoted by ft ra ns ~ 3 the 
same-side 2D peak amplitude and pseudorapidity width 
cr^ A and away-side ID peak amplitude transition to a 
qualitatively new centrality trend. The slopes of the 
trends on centrality measure v increase within one cen- 
trality bin by factors 3.5 to 5. The transitions for 62 and 
200 GeV Au-Au data are located at similar centralities 
as measured by mean participant path length v employed 
as an A-A geometry parameter. The underlying mecha- 
nism for the transition and its relation to other collision 
parameters is presently under study. Within the minijet 
context the large increase in jet-like correlations in more- 
central Au-Au collisions is inconsistent with expectations 
of strong minijet quenching as described in relativistic 
Boltzmann transport and diffusion theories. 

Angular correlations from the present analysis com- 
bined with pQCD predictions of mean jet number per 
A-A collision were used to estimate jet fragment yields 
which agree quantitatively with pQCD calculations of 
parton fragment yields and with measured spectrum 
hard-component yields, implying that about one third 
of the hadronic final state in 200 GeV central Au-Au 
collisions is contained in resolved jets. 

A scaling analysis of event-wise mean-p t fluctua- 
tions revealed that the same-side 2D peak structure in 
p t (rather than number) angular correlations increases 
smoothly with centrality (i.e. no sharp change in slope) 
from below t'trans to well above it. In the minijet context 
the pt angular correlations trend suggests that the SS 
2D number-correlation peak, although strongly modified 
above the transition, is still initiated by semihard parton 
scattering which increases as least as fast as the binary- 
collision trend. A successful theoretical description of the 
correlation structures presented here must not only ac- 
count for the rapid slope and magnitude changes in the 
SS 2D peak amplitude and t;a width, but also the rel- 
atively smooth evolution of the 2D peak azimuth width 
and same-side pt angular correlations above and below 

"trans • 

The A-A collision-energy dependence of the SS 2D 
peak amplitude from the present number-correlation 
analysis is consistent with that inferred for the corre- 
sponding structure in p t angular correlations and with 
the energy dependence of the nonjet azimuth quadrupole 
(as measured by equivalent statistical quantities). All 
follow the log(y / iNN _ ) trend that might be expected for a 
QCD scattering process. 

Two-dimensional histograms were examined for inde- 
pendent "higher harmonics" (multipoles) represented by 
v m with to > 3. Relative to the fit model of Eq. ^ no 
such structure was observed to the statistical limits of the 
data. The same-side 2D peak can be decomposed into a 
Fourier series of azimuth multipoles with a common large 
curvature on pseudorapidity. v m values inferred from ID 
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analysis of projected 2D angular correlations can there- 55 
fore be predicted from the fitted same-side 2D Gaussian 
parameters. Fits including an additional sextupole model 
element for the more-central collision data were shown to 
be equivalent to an alternative model of the same-side 2D 
peak consisting of a 2D Gaussian (with modified param- w 
eters) plus an additional ID Gaussian on azimuth. The 
sum of those two model elements is systematically equiv- 
alent to the fitted single same-side 2D Gaussian. Ad- 
dition of higher multipoles to the 2D model of Eq. Q 
results in model redundancy and ambiguous fit results. 

Results from the present analysis were compared with 
the expected trends from three scenarios for nuclear col- 
lisions at RHIC : (i) Glauber linear superposition of N- 65 
N collisions (the A-A baseline reference), (ii) a locally- 
equilibrated opaque medium with zero mean free path, 
and (iii) an intermediate scenario where partons and 
their hadron fragments randomly scatter in a medium. 
The present analysis reveals substantial deviations from 70 
(i) . The large same-side peak amplitude increase and az- 
imuth width decrease with centrality apparent in Figs. [1] 
and [31 coupled with the smoothly increasing mean-p t 
angular correlation amplitude for v > I'trans and corre- 
sponding p t spectrum hard-component centrality evolu- 75 
tion above the transition, represent the essential experi- 
mental constraints summarized here. 

We are unable to reconcile an interpretation in which 
RHIC nuclear collisions are dominated by minijet struc- 
ture and the observations reported here are described 8 o 
by pQCD with modified fragmentation, with that of a 
strongly-absorptive, even opaque (zero viscosity), colli- 
sion system as described in scenarios (ii) and (iii). If 
the collision system turns out to be effectively opaque 
to few-GeV partons the present observations would be 85 
inconsistent with the minijet picture presented here. 
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Appendix A: Minijets 

Minijets have an experimental and theoretical history 
over more than twenty years (26l. [27ll85j . A growing num-105 
ber of experimental results indicate that minijets make 



significant contributions to the transverse dynamics of 
nuclear collisions above ^/snn ~ 13 GeV [U [H, [36|, 
HI, H3, [j3|. However, the phenomenological definition 
of minijets and their interpretation as true pQCD jets 
is not widely recognized within the heavy ion commu- 
nity. We review the experimental properties of minijets, 
their relation to QCD theory and their manifestations in 
combinatoric correlation analysis as in the present study. 



1. Minijets and calorimeter experiments 

The minijet concept emerged experimentally at the 
Spp S from a UA1 analysis of E t structure with an event- 
wise cone jet finder down to exceptionally small inte- 
grated E t (5 GeV) HH. The resulting E t "clusters" were 
compared with pQCD predictions to test how low in en- 
ergy apQCD jet description is applicable to E t struc- 
ture da [23. 

The UAl analysis determined that E t clusters follow 
an approximate pQCD power-law parton spectrum down 
to 5 GeV. Azimuth correlations between clusters exhibit a 
peak at ir radians expected for back-to-back parton scat- 
tering. The 5 GeV cutoff in calorimeter data was related 
to a 3-4 GeV parton energy equivalent, the difference con- 
tributed by the underlying event [26|. UAl concluded: 
"...one can usefully define jets down to i?t, m i n ~ 5 GeV 
[3-4 GeV parton energy]. ... The agreement of the in- 
clusive [minijet] cross section with [p]QCD over several 
orders of magnitude is quite remarkable" [HJ . 

More recently, the UAl calorimeter analysis was re- 
peated by the STAR collaboration, confirming a jet spec- 
trum consistent with NLO QCD predictions for event- 
wise reconstructed jets down to 5 GeV energy (back- 
ground corrected to 3-4 GeV) [55]. The same pQCD 
parton spectrum was used to describe spectrum hard 
components as fragment distributions in p-p and Au-Au 
collisions [25| . 



2. Minijets and theory 

In Ref. [27j minijet production was considered in a 
pQCD context for anticipated RHIC U-U collisions based 
on UAl observations. "The observed [UAl minijet] rate 
is in agreement with [p]QCD and is quite large." Appli- 
cability of pQCD to minijets (low-p t jets, E t ~ few GeV) 
was studied in detail down to parton jjf" u = 3 GeV/c in 
Ref. [26|. "...a theoretical cutoff of p" im ~ 3 GeV seems 
to describe the observed total minijet cross section with 
E^\Ef w ) > 5 GeV [3 GeV parton energy]." Minijets 
and cross sections in p-p and p-p collisions were also con- 
sidered in [86]. 

The minijet-based Monte Carlo hijing was developed 
specifically to study the role of minijets in p-p and A-A 
collisions. The parton spectrum is given a lower cutoff pq 
with default value 2 GeV/c. hiji ng p redictions are com- 
pared to p-p collision data in (20l . l87l ] where they quanti- 
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tatively match measured minijet correlations (same-side 
amplitude, widths, away-side ridge) [37], HI]. HIJING 
with "jet quenching" disabled applied to A-A collisions 
is equivalent to a Glauber linear-superposition reference 
[Eq. (0] as shown by the dash-dotted curve in the left- eo 
most panel of Fig. |9] hijing with hard parton scattering 
disabled shows no minijet correlations. 

In a recent study of pQCD applied to Au-Au colli- 
sions, measured pt spectrum hard components identified 
with minijets were described quantitatively by pQCD 65 
calculations [25j . The pQCD predictions are compat- 
ible with measured minijet transverse rapidity correla- 
tions on (yti,yt2) [11], especially a parton fragment spec- 
trum mode at p t = 1 GeV/c in both spectra and correla- 
tions. The direct comparison between pQCD and spec- 
trum data confirms a 3 GeV/c parton spectrum cutoff. 

Theoretical descriptions widely assume that minijets 70 
are rapidly thermalized in RHIC collisions, contribut- 
ing to QGP formation. For example, "minijets... will 
be reprocessed by the system and not emerge from 
it" (27[. Thermalization by transport processes is stud- 
ied in [3l], [Hj]. Thermalization time is estimated as 4- 
5 fm/c with T ~ 200 MeV [H|. Within those mod- 75 
els minijet angular correlations should be strongly sup- 
pressed in amplitude, broadened on both 77 and cf>, and the 
scattered parton p t strongly dissipated among final-state 
hadrons. The correlation structures reported here and in 
Refs. [2ll. I36I4381 l43l . l74j . if generated by semihard parton 80 
scattering and fragmentation, contradict the assumption 
of minijet thermalization. 



3. Minijets and combinatoric analysis 



from jet-like angular correlations [62|. That analysis 
combined pQCD estimates of jet number in Au-Au col- 
lisions with factorization of jet-correlated pair numbers 
to obtain parton fragment yields. The results are in ex- 
cellent agreement with fragment yields inferred from a 
two-component analysis of identified-hadron spectra [64| . 

A combination of several related analyses has estab- 
lished a quantitative connection among spectrum hard 
components, pQCD-predicted fragment distributions and 
jet-like correlations, providing strong support for a mini- 
jet interpretation. The combined analysis [o2| demon- 
strates that one third of the hadronic final state in cen- 
tral 200 GeV Au-Au collisions lies within resolved jet 
correlations. 



Appendix B: 2D Fit-model ambiguities 

In Sec. IVIII HI we emphasized multipoles inferred by 
projecting the entire SS 2D peak onto ID azimuth and 
evaluating the SS peak Fourier components. Those mul- 
tipole amplitudes may correspond to "nonflow" bias in 
v m data inferred from nongraphical numerical meth- 
ods (79|. In this Appendix we address a more subtle 
problem: what is the consequence if a sextupole (113) or 
higher multipole element is added to the standard 2D 
fit model used for the present analysis. Is the extended 
model superior to the original? Is more information ex- 
tracted from the data? To answer such questions we 
must distinguish between necessary and sufficient condi- 
tions, lower bounds vs upper bounds on fit parameters 
and algebraic equivalence between seemingly different fit 
models. 



Minijet structure has also been observed directly in 
minimum-bias two-particle angular correlations from 200 
GeV p-p collisions [37], HH . Angular correlations with no 
"jet" minimum-p t requirement exhibit just the structure 
expected from pQCD jets: a narrow intrajet same-side 2D 
peak at the angular origin (parton fragmentation) with 
most-probable pt ~ 1 GeV/c and an interjet away- side go 
ridge at n radians (back-to-back parton scattering). 

The intrajet 2D peak is interpreted to reflect the angu- 
lar consequences of a parton fragment momentum distri- 
bution along the jet axis with mode near 1 GeV/c [25l[75j 
and a soft momentum spectrum transverse to the jet 95 
axis [89j . The resulting 2D "jet cone" typically falls 
within one radian, with single-particle r.m.s. width ~ 0.5. 
The interjet ID peak on azimuth is uniform over a 
large interval on r]A because the parton-parton CM is 
broadly distributed relative to the nucleon-nucleon CM. 100 
The away-side azimuth width is determined by a combi- 
nation of the intrajet width ~ 0.5 and acoplanarity due to 
parton intrinsic k t in the projectile nucleons [37l.[3a]. The 
observed systematics and angular correlations are quali- 
tatively consistent with the PYTHIA Monte Carlo 90]. 105 

Preliminary results from the present analysis [43j were 
used to derive absolute spectrum hard-component yields 



1. Jet-related vs nonjet structure 

We first establish a distinction between "nonjet" and 
"jet-related" structure in 2D angular correlations. Based 
on correlation systematics below the sharp transition at 
vtrans we interpret two elements of the data model in 
terms of (mini)jets: the SS 2D peak and the AS dipole. 
Whether a jet interpretation for those structures is appro- 
priate in more-central Au-Au collisions is an open ques- 
tion. To provide consistent terminology for this analysis 
we refer to those two model elements as "jet-related" 
structure for all collision centralities. All complementary 
structure is by definition "nonjet," including the soft- 
component ID Gaussian on 77A and the BEC-electron- 
pair peak, but also the nonjet quadrupole component 
Aq which is observed to have negligible curvature on t]a 
within the TPC angular acceptance. 

In describing multipole elements v m we distinguish be- 
tween (a) v m obtained from ID Fourier fits to all angu- 
lar correlations projected onto ID 0a and (b) v m ob- 
tained from model fits to 2D angular correlations which 
may include one or more multipole terms Ax (letters 
X = D, Q, S, O denote azmuth multipoles: dipole, 
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quadrupole, sextupole, octupole, with pole number 2m). 
Multipoles in case (a) must mix jet-related and non- 
jet structure as discussed in Sec. IVIII HI Multipoles in 55 
case (b) may represent jet- related and/or nonjct contri- 
butions, as discussed in this Appendix. 

Given those definitions and no significant sextupole 
amplitude in fit residuals from our standard 2D model 
fits we conclude no detectable nonjet (^A-uniform) v% eo 
structure in the data. A test for "false negatives" can 
be performed by adding a known sextupole component 
to the data and refitting with the standard 2D model. 
We then observe a sextupole component in the residuals 
with about 50% of the input amplitude. We therefore 6 5 
conclude from the data that the upper limit on any non- 
jet sextupole is about 10% of the SS 2D peak sextupole 
component, which is negligible compared to other corre- 
lation structure. Thus, the only source for any inferred 
sextupole component is the jet-related SS 2D peak, as 70 
noted in Sec. IVIII HI We now consider the significance 
of possible jet-related sextupole amplitudes in 2D data 
histograms inferred by adding a cos(3</>a) term to the 
standard 2D data model. 



2. SS 2D peak models vs peak properties 

It is important to distinguish between model- 
independent peak properties and peak-related model pa- 
rameters. Model parameters should be accurately con- 
strained by peak properties. But parameters from some 
models may be poorly constrained by data obtained 
within a limited acceptance. In more-peripheral collisions 
the entire SS 2D peak is resolved within the STAR TPC 
angular acceptance. The peak model is then unambigu- 
ous: a 2D Gaussian. The SS peak continues to be well 
resolved until well above the sharp transition in peak 
properties at v tTans . However, in more-central collisions 
the SS 2D peak extends sufficiently far outside the TPC r\ 
acceptance that its modeling on ryA may become ambigu- 
ous. When a fraction of the peak is observed some peak 
properties (e.g., amplitude and curvature at the mode) 
are more reliably determined than others (e.g., higher 
moments, tail structure). Model parameters relating to 
the latter may suffer from large systematic uncertainties. 80 

The curvature of a peak near its mode is a well-defined 
algebraic quantity obtained directly from data as the 
second derivative (second difference) of histogram data. 
The curvature of a Gaussian function Aexp(— x 2 /2a 2 ) is 
A/a 2 . For a single Gaussian fitted to a single ID peak the 
Gaussian amplitude and width parameters are directly 
and accurately related to the peak amplitude and cur- 
vature at the mode, the two most well-determined peak 
properties in all cases. Introduction of additional model 
parameters may result in large systematic uncertainties *> 
if corresponding peak properties (e.g., higher-order mo- 
ments) are not well-determined by the data. 



3. Model equivalence: sextupole vs ID Gaussian 

The model elements Ad (dipole) and Aq (quadrupole) 
already present in the standard 2D data model are in ef- 
fect parts of a (truncated) Fourier series. Adding an As 
(sextupole) term to the data model extends the truncated 
Fourier series. The next term Ao (octupole) is typically 
at the level of statistical fluctuations [73, Ho| ■ The sep- 
arate sinusoids may serve both as representatives of dis- 
tinct (nonjet) data multipoles and as parts of a Fourier 
series representing a localized (possibly jet-related peak) 
structure on azimuth [59l [79| . 

If a sextupole component uniform on tja is added to 
the standard 2D fit model nonzero sextupole amplitudes 
are indeed inferred from model fits to more-central Au- 
Au data, and some other model parameters are shifted 
substantially from their original values. We demonstrate 
below that addition of the As term is equivalent to mod- 
ifying the SS 2D peak model by adding a SS ID Gaussian 
narrow on </>a and uniform on r?A- An illustration is pro- 
vided in Table U for 9-18% central 200 GeV data. The 
argument follows those presented in Refs. [5^, [7^| and in 
text relating to Fig. [5] that equate a truncated Fourier 
series to a periodic peak array. Other mid-central data 
give similar results. Fit parameters that do not change 
significantly and are not relevant to this discussion are 
omitted. The parameter labels are as defined in Eq. (QJ. 

TABLE I: Model parameters for three fit models: (a) standard 
model for this analysis, (b) standard model plus additional 
sextupole term As and (c) standard model plus additional 
SS ID Gaussian on 0a with amplitude Aid as part of the SS 
2D peak model. The fit parameters are as defined in Eq. Q. 
Uncertainties in the second column illustrate typical fit uncer- 
tainties (statistical plus systematic) for each parameter. Ax 
denotes the additional model parameter As or Aid- 



parameter 


standard 


std + A s 


std + Aid 


A 1 


0.76 ± 0.04 


0.51 


0.47 




2.3±0.3 


1.78 


1.72 


Aq 


0.18±0.008 


0.23 


0.18 


Ad 


0.29±0.02 


0.175 


0.28 


A x 




0.014 


0.29 



The second column (standard) shows fit results ob- 
tained with the standard 2D model. The third column 
(std + As) shows fit results when a sextupole term 
As cos(34>a) is added to the 2D model and all else re- 
mains the same. Proof that addition of a sextupole term 
to the standard 2D model is equivalent to adding a ID 
Gaussian on azimuth to the SS 2D peak model proceeds 
as follows. 

Table |TT] (third column, data) shows the changes in 
fitted multipole amplitudes between the standard model 
and that including the sextupole term. The third column 
(F m defined in Sec. IVIII Hp shows the calculated Fourier 
coefficients for a unit-amplitude ID Gaussian on azimuth 
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TABLE II: Comparison of changes in multipole amplitudes 
when an additional sextupole element is added to the stan- 
dard model with Fourier coefficients of a SS ID Gaussian on 
azimuth. The comparison demonstrates an equivalence be- 
tween the changes in multipole amplitudes (columns 2 and 3 
of Table U} and the Fourier coefficients of a ID Gaussian on 
azimuth (column 5 of this Table). Column 5 is obtained by 
multiplying Fourier coefficients F m by the value Aid = 0.25 
that best matches the third column. That value also corre- 
sponds to the difference AAi = 0.76 - 0.51 from Table U 



m 


parameter 


data 


i m 


F m Am 


1 


-AA D 


0.115 


0.44 


0.11 


2 


AA Q 


0.05 


0.21 


0.052 


3 


AA S 


0.014 


0.062 


0.0155 



factorization accurately describes the SS peak in all cases, 
and we can simplify the model choices on 7/a to a single 
ID Gaussian with parameters (A, o~ n ) vs a ID Gaussian 
plus constant offset with parameters (A', a' Aid). 



0.15 = 




0.05 - 



with width =0.7 [7^|. The width value was adjusted 
to achieve the best match between columns 3 and 5. The 
fifth column {F m Aio) shows the Fourier coefficients for 
a SS ID Gaussian with amplitude Aid — 0.25. Compar- 
ison of the third and fifth columns reveals equivalence 
within small data uncertainties. The inferred amplitude 
Aid of the added SS ID Gaussian on azimuth corre- 
sponds to the difference AAi between SS 2D peak am- 
plitudes for the two fit models. And the inferred width 
corresponds to the azimuth width of the SS 2D peak in 
the standard fit model. Thus, the parameter changes re- 
sulting from inclusion of a sextupole element in the 2D 
model are equivalent to adding a constant offset to the 
7/a factor of the SS 2D peak model. 

Table U (fourth column) confirms that equivalence by 
presenting the result of a 2D fit replacing the added sex- 
tupole element by an offset in the t]a factor of the SS 2D 
peak model [A x exp(-r&/2tf A ) -> A' x cxp(-n 2 A /2a' 1 f A ) + 
Aid]- The added constant in the ?7a factor then rep- 
resents a SS ID Gaussian on 0a in the 2D model [see 
Eq. (0}]. The fit results are systematically consistent 
with the previous exercise: the ID Gaussian amplitude 
is Aid = 0.29 compared to 0.25, and the modified 
2D amplitude A[ = 0.47 is consistent with Aid (i-e., 
A[ + Aid = A\ = 0.76). Modification of rj widths in the 
revised models is discussed below. Note that Ad and Aq 
have reverted to values for the standard fit model, since 
the added ID Gaussian on azimuth is modeled explicitly 
as a Gaussian rather than as a truncated Fourier series. 
Equivalence of a ID periodic peak array and a truncated 
Fourier series on azimuth is discussed in Refs. I59L [731. 



4. Model ambiguities and parameter significance 

We now return to the issue of SS 2D peak modeling 
and the significance of model-fit results. The standard 
model of the SS 2D peak for this analysis in Eq. (|4]) 70 
is factorized to ID functions on 7/ a and 0a- Whatever 
the actual peak structure on 77A we do observe that the 
narrow ID Gaussian on 0a is independent of ?7a- Thus, 



FIG. 10: Left: The data histogram is the 2D histogram in 
Fig. (right panel) projected onto ID t/a- The narrow peak 
at the center is BEC + electron pairs. The dash-dotted curve 
through the data histogram is a single ID Gaussian with pa- 
rameters (A, a v ). The dashed curve is the sum of a ID Gaus- 
sian (A', a' v ) plus constant offset Aid- The two curves are 
barely resolved at the r\ accepance limits. The dotted line is 
the constant offset. Right: The second difference of the data 
histogram in the left panel corresponds to its local curvature. 
The calculated curvature A/ai of the ID Gaussian in the left 
panel (dash-dotted curve) is the dash-dotted line in this panel. 

Figure [10] (left panel) shows the histogram in Fig. [7] 
(right panel) projected onto ID r/A- The projection over 
27r azimuth eliminates the small tja modulation of the 
AS dipole discussed in Sec. I VII CI Those data can be de- 
scribed in a model-independent way with a Taylor series. 
Two series terms are systematically significant within the 
77 acceptance: the amplitude (constant) and the curva- 
ture (quadratic). The quartic term is not significant. 
The data are thus accurately described by a ID Gaussian 
(dash-dotted curve) with amplitude A — 0.195 and width 
(j f; = 2.9. The Gaussian curvature is Aj a 2 = 0.023. 
In the right panel the Gaussian curvature value (dash- 
dotted line) is plotted together with the second difference 
of the data histogram, comparing local statistical noise 
to the curvature signal. The fourth difference (quartic) 
would be overwhelmed by statistical fluctuations. 

If a constant offset is added to the ID fit model the data 
can constrain only certain parameter combinations. The 
constraints are A = A'+A\d for the common peak ampli- 
tude and A/cr 2 = A' jo' 2 for the common peak curvature 
at the mode. The parameter combination (A 1 , a', Aid) 
is otherwise free to vary within those constraints. If we 
choose the value Aid — 0.077 (dotted line, equivalent to 
Aid = 0.29 for the unprojected 2D peak model) the data 
constraints determine that A' = 0.12 and a' v = 2.27. A 
dashed curve representing the sum of the modified ID 
Gaussian and the constant offset is just visible in the left 
panel above the dash-dotted curve at the r\ acceptance 
limits. If o-' v is reduced by 5% the curves are not visu- 
ally distinguishable. The same exercise could be carried 
out with a large range of Aid choices including zero (the 
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standard data model). Returning to the 2D fit parame- 
ters in Table U we have A = 0.76, A' = 0.47, A 1D = 0.29 
and a v = 2.3. The calculated er^ = 1.82 when reduced 
by 5% is 1.73, consistent with the free- fit value 1.72 (see 
Table U. 

SS 2D peak rj structure in data obtained within a lim- 
ited rj acceptance, as in the present case, can only con- eo 
strain two model parameters. Introduction of the offset 
Aid (i-e., the ID Gaussian on azimuth) as an additional 
model element leads to large fit instabilities because the 
data within the limited rj acceptance cannot constrain 
that parameter. The actual Aid values from model fits 
depend on small statistical fluctuations and systematic 65 
distortions at the 77 acceptance boundary, as in Fig. [10] 
(left panel). The additional model parameter acts to 
amplify relatively small and systematically insignificant 
variations in the data to appear as relatively large but 
still systematically insignificant variations in an extra- 70 
neous model parameter. If the unnecessary element is 
introduced as a Fourier component (t'3) the definitions 
of other multipole amplitudes change substantially, and 
those parameters may then be misinterpreted. 

If a larger fraction of the SS 2D peak is accessible 75 
within a detector angular acceptance more peak prop- 
erties are determined, such as higher peak moments 
(e.g., kurtosis) and possible non-Gaussian tail structure 
at larger rj^- In more-peripheral Au-Au and p-p colli- 
sions essentially the entire SS 2D peak is resolved within 
the STAR TPC acceptance, and the 2D peak model is 
fully constrained. If the detector rj acceptance were ex- m 
tended then even in more-central Au-Au collisions peak 
structure at large ?7a might be accurately determined. 



5. Related systematic uncertainties 
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Appendix IB 31 describes the effects of adding unneces- 
sary model elements to a 2D fit model, and demonstrates 
the equivalence of adding either a W3 or SS ID Gaussian 
element to the model presented in this paper. The nu- 
merical analysis in that subsection demonstrates that any 
changes in the actual data description are less than 3% 
of the reported parameter values. For instance, A\ and 
the sum A' x + Aid agree to better than 3%. Also, when a 90 
SS ID Gaussian element is added to the standard model 
the fitted Ad and Aq amplitudes change by less than 2% 
of their values. 

We further supplemented the systematic uncertainties 
analysis in Sec. IVII El with a consideration of the pos- 95 
sible effects of an undetected U3 component or SS ID 
Gaussian in the data. Data from representative central- 
ity bins 38-46% and 9-18% for 200 GeV collisions were 
studied. 2D fit residuals from the standard model were 
fitted with either a V3 or SS ID Gaussian model com-100 
ponent. Inferred amplitudes As were less than 0.0003 
and inferred amplitudes Aid were less than 0.01. Cor- 
responding additional model elements with the inferred 
amplitudes held fixed were subtracted from the standard 



data model and the data were refitted with the revised 
model. Changes in the standard-model parameters were 
in all cases substantially less than the statistical fitting 
errors reported in this paper. 



6. Summary 

Ideally, data models should be composed of elements 
that are both necessary and sufficient. Removal of a nec- 
essary element from a data model results by definition 
in appearance of equivalent structure in the fit residu- 
als. The model is incorrect without that element. A 
combination of necessary elements that leaves no signifi- 
cant structure in the fit residuals is sufficient. A model 
composed of a sufficient assembly of necessary elements 
accurately and efficiently describes the data. The model 
defined by Eq. ((4]) is a necessary and sufficient data model 
for 2D histograms obtained within the STAR TPC an- 
gular acceptance, except as noted in Sec. IVII CI 

Within the STAR TPC acceptance we find that offset 
Aid values inferred by adding a SS ID Gaussian on az- 
imuth to the standard 2D fit model, or equivalently U3 
values inferred by adding a sextupole element, are not 
systematically significant. The results of this analysis 
have no implications for SS 2D peak structure at larger 
?7a outside the STAR TPC angular acceptance. 



Appendix C: Constructing f a b 

Pair ratio r = p S ib/prcf is measured most accurately 

for a single charge combination (i.e., H — h, , or +— ), 

event multiplicity N c h and vertex position z v tx within the 
STAR TPC [iU because in that case charge-dependent 
tracking inefficiencies and acceptance effects most closely 
cancel in the ratio. In a practical analysis particle pairs 
from small bins on iV c h (subdivisions of a centrality bin) 
and z v tx are combined into individual ratios denoted by 

f'a.ab — ^a,ab,sih/ ^a,ab,rcf ) (^-0 

where subscript a represents bin indices for charge-pair 
combination, event-multiplicity and vertex-position bins. 
Indices (a, b) denote (?7a,<M) bins. We define normal- 
ized pair numbers n a ^ ab = n a ^ ab / J2 a b n a,ab for sibling or 
mixed (ref) pairs, where n a ,afc is the ensemble-averaged 
event- wise number of pairs in 2D angle bin (a, b), with 
additional bin indices a. We obtain stable (optimized) 
angular correlations when 7V c h and z v tx bins are restricted 
to bin widths AiV c h < 50 and Az vtx < 5 cm respectively. 

Corrections were applied to each ratio r a ,ab for two- 
particle reconstruction inefficiencies due to overlapping 
space points in the TPC (two trajectories merged into one 
reconstructed track) and intersecting trajectories which 
cross paths within the TPC and are reconstructed as 
more than two tracks (splitting) (50j . The track- merging 
cuts rejected track pairs if both the longitudinal (along 
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the TPC drift direction) and transverse (along the pad 
row direction) separation distances in the TPC were 
less than 5 cm at any of three radial positions in the 55 
TPC [42[ (inner, middle and outer radii). The crossing 
cut rejected track pairs if their longitudinal separations 
at the middle and outer radii in the TPC were both less 
than 5 cm and if they cross in the transverse plane, i.e. 
(01, inner 02, inner) (01, outer -02, outer) < 0, where angles 60 
for tracks 1 and 2 locate trajectory azimuth positions in 
the TPC at the inner or outer radii. The same track-pair 
cuts were applied to sibling and mixed pairs. Cuts have 
small overall effect except near (t?a,0a) = (0,0) where 
artifacts (depressions in f a b — 1) & re significantly reduced. 

The grand ratio for each angle and centrality bin was 
obtained by constructing a sibling-pair-weighted average 
over f Qja b for the relevant charge combinations (all four 
are summed), multiplicity bins, and vertex-position bins 
indexed by a. The final normalized ratio is defined by 

?ab{v) = 2_j N ot,sib fa iab / 2J ^Va,sib (C2) 65 

a(t>) a (y) 

for centrality bin v, where N a ,^b = J2ab n a,abMb- 

Appendix D: Event pileup corrections 

Event pileup can contribute spurious structure to an- 
gular correlations. Ideally, each triggered collision event 
would appear in isolation in the TPC and be digi- 
tized over a readout time interval of 40 /is (TPC drift 
time [42|). At sufficiently high beam luminosity particle 
tracks from an untriggered A-A collision in a different 
beam-beam crossing (every 120 ns) may appear in the 
TPC tracking volume during readout of a triggered colli- 
sion event, a source of contamination known as "pileup." 
The extraneous particle tracks may be added to the trig- 
gered event or may be misidentified as the triggered col- 
lision event by the reconstruction code. We estimate 
that approximately 0.5% and 0.05% of the minimum-bias 
events in the 62 and 200 GeV Au-Au data sets respec- 
tively are contaminated by pileup. 

Although relatively few in number, pileup events can 
significantly affect angular correlations, particularly the 
r)A dependence. We observe that event pileup produces 
a characteristic "W" correlation shape on tja due to the 
mis-match in r\ range of particle distributions from out- n 
of-time partial events (~1 unit) and those for intact ref- 
erence events (2 units). Due to the ry-range mismatch all 
particles in pileup events appear as correlated pairs rela- 
tive to the normal (two units in 77) reference-pair distri- 
bution, thus greatly amplifying the relative contribution 
of pileup events to the total angular correlations. For ex- 
ample, a 0.5% pileup event rate contributes a few parts 
permil to sibling-to-mixed pair ratio r, which is com- ao 
parable to the true correlation amplitude. The pileup 
structure has no significant 0a dependence. 

Pileup events were identified and removed (with es- 
timated 75% efficiency) by exploiting the bi-directional 



drift of the STAR TPC which causes reconstructed par- 
ticle tracks from pre- and post-triggered collisions to be 
either split at the TPC high- voltage central membrane or 
truncated before the TPC readout plane. These tracking 
artifacts produce distinctive patterns in the event-wise 
track end-point distributions in the longitudinal drift di- 
rection for trajectories within the active volume of the 
TPC, leading to an efficient pileup filter algorithm. 

Pileup-free correlations Ap/^/p le { were constructed 
from 2D histograms obtained with and without the pileup 
filter by solving the equations 

_J^ (nonlter ) = -% + ^(pileup) 

VPrrf y/Prcf y/Pref 

(with filter) = + (1 - f)-^L (pileup), 

(Dl) 

where Ap/ ^Jp rc f (pileup) represents the pileup contribu- 
tion. The centrality-independent pileup detection effi- 
ciency (1 — /) = 0.25 ± 0.1 assumed for this analysis 
is based on an estimate of the fraction of pileup events 
which had too few tracks crossing the central membrane 
to be identified by the adopted pileup signature. The es- 
timated uncertainty in (I — /) is propagated to the total 
uncertainties for the analysis results. 




FIG. 11: (Color online) Left: Uncorrected angular correla- 
tions from 62 GeV 37-46% central Au-Au collisions showing 
pileup distortions, especially evident as the W-shaped nonuni- 
formity of the away-side ridge on r/A- Right: The same data 
with pileup correction applied. 

Figure [TT1 shows uncorrected (left panel) and corrected 
(right panel) histograms for mid-central 62 GeV Au-Au 
data where pileup distortions are most severe. Although 
pileup may affect a small fraction of the total event num- 
ber (e.g. < 1%), its effect on angular correlations can be 
substantial, as illustrated in the figure. 

Appendix E: 130 GeV Au-Au data 

The 2D angular autocorrelation method used in this 
paper was first applied to unidentified charged hadrons 
from Au-Au collisions at 130 GeV [21] . Similar correla- 
tion structures were observed in the 130 GeV analysis, 
but with poorer statistics and coarser centrality binning. 
Here we compare 130 GeV data to the present analysis. 
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The 130 GeV centrality bins were denned by approxi- 
mate cross-section fractions 40-70%, 17-40%, 5-17% and 
0-5%. After conversion to measure Ap/^/p ro f using 
log y'SNN interpolated prefactor d 2 N c h/drjd(/) the same- « 
side 2D Gaussian peak amplitude = 0.14, 0.22, 0.26, 
0.24; a nA = 0.58, 1.05, 1.34, 1.36; ct 0a = 0.61, 0.55, 0.54, 
0.53; dipole amplitude = 0.023, 0.052, 0.067, 0.057; and 
quadrupole amplitude = 0.07^ 0.186, 0.090, 0.025 for fit 
parameters in Table I of Ref. [21| for the four respective 50 
centralities. The quadrupole amplitude at 130 GeV inter- 
polates the present results at 62 and 200 GeV. However, 
the dipole amplitude and the SS 2D Gaussian parame- 
ters do not, each being smaller than expected from linear 
interpolation of the 62 and 200 GeV parameters. 55 

The 130 GeV correlation analysis differed from the 
present analysis in several respects: (a) In Ref. [2l[ a 
Pt upper limit of 2.0 GeV/c was imposed, whereas in this 
analysis there is no imposed upper limit, (b) Pair cuts 
were applied to minimize quantum correlations (HBT) 60 
and conversion-electron contamination near the angu- 
lar origin. Those cuts reduced the correlation ampli- 
tude near (v7a,</>a) = (0,0) by about 20%, including the 
same-side 2D Gaussian peak. In the present analysis no 
such cuts were applied — the HBT/electron contribution 
is explicitly accounted for in the model function (the A? 
term) . (c) The 2D exponential term in the present fitting 
model was not included in Ref. [2l|. (d) TPC two-track 
merging near the angular origin extended to twice the 
angular range of more recent data because the 130 GeV 
data were obtained with half the nominal magnetic field 
strength of the STAR detector, (e) The 130 GeV analy- 
sis accepted a larger pseudorapidity range (|?7| < 1.3) and 
included collision vertices distributed over a much larger 
distance interval (|z v tx| < 75 cm). 70 

It was not feasible to repeat the 130 GeV analysis start- 
ing from the raw data. However, the above p t < 2 GeV/c 
restriction and extended pair cuts were applied to the 
present 200 GeV data in order to estimate their effects on 
the fitting results. The additional p t and pair cuts caused 75 
the 2D Gaussian amplitude and widths to decrease as did 
omitting the 2D exponential term in the fitting model. 



The dipole amplitude was similarly reduced, however the 
quadrupole amplitude was weakly affected except for the 
most-central 0-5% bin where it increased. Omission of 
the 2D exponential had greater effect on the 2D Gaus- 
sian widths than on the amplitudes. The additional pseu- 
dorapidity range, increased track merging and extended 
collision vertex range in the 130 GeV data had negligible 
effects on the fit parameters. 

Reversing the above fractional changes in amplitudes 
and shifts in widths and applying those changes to the 
130 GeV fit parameters results in better agreement with 
expected energy-dependent trends, where: A\ — 0.14, 
0.32, 0.43, 0.31; a, 1A = 1.5, 1.8, 2.3, 1.9; cr 0A = 0.81, 
0.67, 0.65, 0.58; A D = 0.07, 0.11, 0.14, 0.08; and A Q 
= 0.07, 0.17, 0.05, 0.01. The fractional amplitude and 
width increases introduce additional ±10% uncertainties 
beyond that in Ref. [Hj]. 

This estimation procedure is not intended as a substi- 
tute for a re-analysis of the 130 GeV raw data or even 
new measurements at 130 GeV. It does however provide 
a reasonable explanation for the discrepancy between the 
130 GeV parameters in Ref. [2l[ and what would be ex- 
pected given the present results at 62 and 200 GeV. 



Appendix F: Tabulated data 

Model fit parameters, statistical errors (in paren- 
theses), and asymmetric systematic uncertainties (sub- 
scripts and superscripts) for 200 GeV and 62 GeV Au-Au 
charged-particle angular correlation data for eleven cen- 
trality bins from most-peripheral (left column) to most- 
central (right column) are listed in Tables IIIII and IIVI 
The volume of the same-side 2D Gaussian peak within 
the acceptance, corrected multiplicities, and Monte Carlo 
Glauber centrality measures v and N vax t are also listed. 
The full volume of the same-side 2D Gaussian peak ex- 
trapolated to 47T acceptance is given by 27rAi<7 J?A er (; i A . 
X 2 /DoF values are for 158 degrees of freedom. 
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stated limits. angular correlations. 

"Jet-like angular correlations" refers to a generic same- 
side 2D peak and away-side ID ridge structure in 2D 



TABLE III: Model fit parameters, statistical errors (in parentheses), and asymmetric systematic uncertainties (subscripts and superscripts) for 200 GeV Au-Au charged-particle 
angular correlation data for eleven centrality bins from most-peripheral (left column) to most-central (right column). Prefactor (d 2 N c h / 'dr)d<j)) uncertainty (Sec. IVII Aj) is not 
included. The volume of the same-side 2D Gaussian peak within the acceptance, corrected multiplicities, and Monte Carlo Glauber centrality measures v and N paTt are also listed. 
\ 2 values are for 158 degrees of freedom. 
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0.224(21) + '°" 

V / — .017 


Aq 


0.002(l)+°,°,° 5 

^ WVJ ^V / — .001 


0.011(1)1; 


0005 
001 


0.028(1)1 


.001 
.002 


0.070(l)+nni 


0.136(2)+°,°4 

^ V / — .007 


0.201(4)1;°™ 


0.270(4)1 


.006 
.012 


0.268(5)1 


.007 
.012 


0.179(6)1 


.005 
.010 


0.063(9)1 


.001 
.010 


0.001(8)+-°,°,? 

J v w V ./ — .007 


Ai 


0.091(4)+°°! 

V / — .00u 


0.105(10)1 


-.016 
-.015 


0.142(6)1 


.012 
.011 


0.169(6)+S? 

^ "^V"/ — .Ola 


0.207(7)+'n 22 

1 V 1 / — .011 


0.316(19)1 


-.029 
-.020 


0.378(17)l;°g 3 9 


0.590(28)1 


-.039 
-.021 


0.767(38) H 


-.048 
-.031 


0.760(50); 


K054 
-.006 


0.646(40)+'^ 

V / — .012 


'/A 


0.585(40)+°,o 3 

\ / — .031 


0.699(45)1 


-.054 
-.002 


0.748(34)1;°^ 


0.771(35)1 131 


0.892(45)1 355 


1.372(78)1 


-.364 
- 032 


1.470(64)1 


-.562 
-0 


1.972(74)1 566 


2.320(85)1 


- . 589 
-0 


2.185(107) 


+ .561 
-0 


2.134(108)1 656 




0.750(41)1;°^ 


0.974(72)1 


-.050 
-.049 


0.818(28)j 


-.031 
-.030 


0.756(23)1;°! 


0.662(14)1;°^ 


U.6(0( v 16J_ 


-.016 
-.032 


0.626(13)1 


-.009 
-.029 


0.659(11)1 


-.015 
-.010 


0.677(11)1 


-.015 
- .007 


0.663(15); 


K015 
-.001 


0.629(15)1;°°? 


A 2 


0.184(43)+°,°J! 

V / — .009 


0.198(10)1 


-.008 
-.008 


0.335(13);! 


-.009 
-.009 


0.466(19)+°1? 

V / — .011 


0.656(21)+°" 

V " } —.013 


0.848(24)1 


-.018 
-.014 


1.120(33)1 ° 2 * 


1.541(48)1 


-.024 
-.012 


z.loU^y i ) _ 


-.031 
-.011 


4.3±1.4l 


.036 
.008 


3.8(8)+-°« 

V / - .007 




0.059(21)^.054 


0.218(33)1 


-.012 
-.012 


0.164(12)1 


-.001 
-.001 


r\ i 4-i /i a\ + .008 

0.141(10)1^ 


r\ 1 -i o /K\ .008 

0.113(5)+^ 


0.106(4)1 


.004 
.005 


0.090(3)1 


.004 
.004 


0.081(2)1 


.002 
.002 


0.069(2)1 


.001 
.002 


0.058(2)1° 003 


0.056(2)1° 00 4 


W 4>A 


0.126(27)1;°,!* 


0.296(45)1 


-.009 
-.009 


0.183(14)1 


-.008 
-.007 


0.157(10)1;°° 7 7 


0.133(7)1;°^ 


0.119(5)1 


.012 
.011 


0.093(4)1 


.011 
.008 


0.074(3)1 


.004 
.001 


0.056(3)1 


.006 



0.029(11); 


K002 
-.005 


0.036(7)1;°° 2 


A 


0.034(2)+SS? 

V / — .00o 


0.015(2)1;°° 1 1 


0.013(2)1 


.001 
.001 


0.008(2)+-^ 

V J — .002 


0.004(l)+°° 2 

\ / — .00a 
























(TO 


0.593(44)1;°^ 


0.398(56)1 


-.019 
-.019 


0.309(42)1 


-.027 
-.028 


0.118(72)1;°" 


0.210(66)1; 2 ?? 
























A 3 


-0.023(l)l°° 1 i 


-0.025(2)1 


-.001 
-.002 


-0.025(2)1;°° 2 


-0.025(2)1;°°« 


-0.028(3)l;°f 2 


-0.057(7)1 


.008 
.013 


-0.065(6)1;°^ 


-0.124(10); 


f-.009 
-.019 


-0.178(14); 


f .011 
-.020 


-0.172(18) 


-I-.006 
-.019 


-0.142(14)1;°°? 


volume 


0.251 


0.447 




0.542 




0.613 


0.749 


1.561 




1.806 




3.322 




4.627 




4.427 




3.549 


dN ch /dr] 


5.2 


13.9 




28.8 




52.8 


89 


139 




209 




307 




440 




564 




671 


V200 


1.40 


1.68 




2.00 




2.38 


2.84 


3.33 




3.87 




4.46 




5.08 




5.54 




5.95 


-/Vpart 


4.6 


10.5 




20.5 




36.0 


58.1 


86.4 




124.6 




176.8 




244.4 




304.1 




350.3 


X 2 /DoF 


1.00 


1.00 




0.96 




1.22 


1.71 


1.92 




2.96 




3.11 




3.68 




2.00 




2.57 



to 



TABLE IV: Same as Table Hill except for 62 GeV Au-Au charged-particle angular correlation data. 

62 GeV Au-Au 



V AA \/0) 




7^ 8/1 




OO- / 


OD-DO 




40-00 


37-46 




Zo-O ( 


18-28 




y-io 


K Q 

o-y 


u-o 


A D 


0.028(1)1 U,W 


0.038(1)1 


.(104 

.001 


0.042(1)1;°^ 


0.048(2)1 


.005 

.00'' 


0.056(2) l;^ 


0.074(4)1;^ 


0.121(8)1;^ 


0.176(12)1;^ 


0.197(13)1;^ 


0.160(13)l o u "° 


0.153(15)1 UO " 




0.004(l)±;g£ 


0.007(1)1 


.001 
001 


0.016(1)1;°°? 


0.037(1)1 


.002 
002 


0.073(1)1;°°^ 


0.117(1)1 


.007 
.004 


0.147(2)± : «g 


0.148(4)1 


.008 
.004 


0.106(4)1;° ° 


0.048(4)1 014 


0.003(5)1 ° 16 


Ai 


0.043(3)1;°^ 


0.072(3)1 


.004 
(ION 


0.084(3)1;°°]; 


0.109(3)1 


.007 
Oil 


0.128(4)1;°?^ 


0.176(8)1 


.006 
.025 


0.284(16)1;^ 


0.407(24)1 


-.024 
-.047 


0.458(25)l;os8 


0.389(26)l° 0g4 


0.367(30)l° lo2 




0.559(40)i;°^ 


0.666(28)1 


-.027 
- 018 


0.720(27)1;^ 


0.794(26)1 


-.038 
- 029 


0.907(36)1;°^ 


1.115(52)1 


-.170 
-.128 


1.567(69)1;^ 


1.981(84)1 


-.237 
-.126 


2.104(80)1;^° 


2.120(105)1;^ 


2.159(128)1;|°° 




1.042(82)1;°^ 


1.039(50)1 


-.033 
- 036 


1.001(46)1;°^ 


0.859(22)1 


-.035 
- 040 


0.787(15)1;^ 


0.736(12)1 


-.028 
-.038 


0.727(12)1;°^ 


0.741(12)1 


-.012 
-.026 


0.725(11)1;°°^ 


0.672(15)1° 058 


0.689(18)1° 064 


A 2 


0.086(9)1;°°! 


0.152(6)1 


.007 
.006 


0.251(7)1;°°; 


0.375(9)1 


.008 
.008 


0.495(10)1;°*° 


0.689(13)1 


-.013 
-.011 


0.935(16)1;°^ 


1.223(24)1 


-.020 
-.010 


1.676(40)1;°,?? 


2.360(170)1;°^ 


4.2±1.0l;gg? 




0.100(19)±;^° 


0.175(14)1 


-.022 
-.022 


0.169(9)1;°°« 


0.136(5)1 


.004 
.005 


0.133(4)1;°°^ 


0.113(3)1 


.001 
.002 


0.098(2)l;°°i 


0.086(2)1 




.003 


0.075(1)1° 002 


0.072(2)1° 005 


0.063(2)1° 005 


W 0A 


0.227(35)_ 050 


0.223(18)1 


-.022 
-.022 


0.195(,1UJ_ 012 


0.152(5)1 


.004 
.004 


0.134(4) _ 004 


0.114(3)1 


.006 
.004 




0.075(2)1 


.005 





U.U4z(oJ_Q 


U.Uzo(lUJ_o 


A 


0.036(2)1;°°^ 


0.024(2)1 


.002 
.002 


0.017(1)1;^ 


0.013(1)1 


.001 
.002 


0.013(1)1;°" 


0.007(1)1 


.002 
.002 


- 






- 


- 


- 


cro 


0.625(26)1;^ 


0.480(23)1 


-.018 
-.018 


0.398(31)1;^ 


0.347(30)1 


-.053 
-.057 


0.233(20)1.026 


0.153(25)1 


-.022 
-.029 


















A 3 


-0.021(1)1 °° 2 


-0.021(1)1;°°? 


-0.021(1)1;°°? 


-0.022(1)1;°°? 


-0.023(2)1;° 00 °. 


-0.031(3)l;°i« 


-0.059(6)1;°^ 


-0.097(9)1;°?° 


-0.110(9)1° 2 04 


-0.092(9)1 ° 29 


-0.090(11)1;°^ 


volume 


0.157 


0.312 




0.378 


0.462 




0.558 


0.841 




1.623 


2.580 




2.889 


2.279 


2.215 


dN ch /drj 


3.3 


9.3 




19.8 


36.9 




62.8 


99.0 




149 


217 




312 


403 


488 


f62 


1.34 


1.59 




1.87 


2.22 




2.62 


3.05 




3.49 


3.96 




4.46 


4.82 


5.13 


AT part 


4.0 


9.3 




18.7 


33.8 




55.7 


84.8 




122.0 


171.6 




238.4 


297.6 


344.6 


X 2 /DoF 


1.14 


1.30 




1.08 


1.51 




1.26 


1.65 




2.18 


2.72 




3.56 


3.63 


3.32 



